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ABSTRACT 


' 

An investigation was made of the mobility and types of point defect intro- 
duced in platinum by deformation in liquid nitrogen, quenching into water 
from 1600°c, or reactor irradiation at 50°c. In all cases the activation energy 
for motion of the defect was determined from measurements of electrical 
resistivity. Measurements of density, hardness, and x-ray line broadening 
were also made where applicable. 

These experiments indicated that the principal defects remaining in 
platinum after irradiation were single vacant lattice sites and after quenching 
were pairs of vacant lattice sites. Those present after deformation in liquid 
nitrogen were single vacant lattice sites and another type of defect, perhaps 
interstitial atoms. 


$1. IyTRODUCTION 


ALTHOUGH an equilibrium number of the point defects, vacancies and 
interstitial atoms are present in a metal, a larger concentration is thought 
to be introduced by cold working, quenching, or irradiating. The con- 
centration, distribution and type of the extra defects retained in the sample 
will depend on the temperature and method of production. If the defects 
are mobile at the temperature of production, they will move about the 
lattice until they reach a position of lower energy where they are trapped 
or annihilated. Foreign atoms are expected to act as trapping centres, 
while grain boundaries voids, free surfaces, and perhaps dislocations 
will act as annihilation centres or sinks. An interesting review on the 
effect of lattice defects on some physical properties of metals was recently 
given by Broom and Ham (1958). 

Many different physical properties have been measured to detect 
vacancies and interstitital atoms. Electrical resistivity, which is very 
sensitive to their presence, is often used to measure their mobility. 
From this the activation energy required for one atomic jump of the 
defect can, be calculated and the type of defect can be inferred. Also, the 
number of jumps required by the defect to reach the sink can be deter- 
mined, from which the type of sink can be inferred. Unfortunately, all 


point defects increase the resistivity so that changes in other physical 
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properties are necessary to determine exactly what type of defect is present. 
Changes in density and lattice parameter should discriminate between 
vacancies and interstitial atoms. The x-ray line width is unchanged by 
point defects but indicates the amount of lattice strain produced by 
dislocations in deformed samples. 

In the present experiments the point defects were introduced by cold 
work at —196°c, by quenching from above 1600°c, and by fast neutron 
irradiation at 50°c. The electrical resistivity, density, lattice parameter, 
hardness, and x-ray line width were measured before and after the 
defects were introduced into the sample. In, addition, the mobility of the 
defects was determined from changes of electrical resistivity due to 
heat treatment. 


§ 2. CHOICE OF MATERIAL 


The choice of sample material for these experiments was mainly deter- 
mined by the irradiation conditions. The sample should have a high 
melting point in order to retain a maximum amount of radiation damage 
at 50°C, the temperature of irradiation. In addition it should have a low 
gamma activity after irradiation. These two conditions, combined with 
the usual requirements of simple crystal structure, high purity, and good 
surface properties, made platinum the obvious choice. 

The platinum used was Johnson Matthey and Mallory thermocouple 
grade (99-99%), so that direct comparison could be made with the 
platinum self-diffusion experiments of Kidson and Ross (1957), which 
were made on similar material. 


§ 3. MeTHops For INTRODUCING DEFECTS INTO PLATINUM 
3.1. Deformation 


The 0-010in. diameter platinum wires were annealed in air at 900°C. 
for 30 min, and then extended 8% in liquid nitrogen. They were removed 
from the extension apparatus while immersed in liquid nitrogen. Since 
the potential leads were silver-soldered to the sample before the extension, 
the resistivity could be measured after any recovery anneal above 
— 196°C. 

The hardness samples were 0-lin. thick platinum flats rolled to 8% 
reduction, with the sample immersed in liquid nitrogen between each pass. 


3.2. Quenching 


An electric current was used to heat a 4in. long 0-010 in. diameter 
platinum wire. The potential leads of 0-003 in. diameter platinum wire 
were spot-welded one inch from the ends of the sample. The tempera- 
ture of the heated sample was determined from its resistance. The 
calibration curve for resistance against temperature is shown in fig. 1. 
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The sample was cooled rapidly (quenched) by sudden immersion into 
cold water. To reduce the deformation during this quench, the axis of 
the wire sample was perpendicular to the surface of the water. The time 
taken to cool the sample was measured by retaining the heating current 
and recording the voltage across the sample on an oscillograph. A typical 
quench, shown in fig. 2 (Pl. 25), gives a time of 0-02sec to cool from 
1600°c to 20°c. 


Fig. 1 
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Electrical resistivity of thermocouple grade platinum. No correction was made 
for the change in sample dimensions with temperature. 


3.3. Fast Neutron Radiation 


The 0-010 in. diameter platinum wire was annealed in air at 900°C 
for 30min and then irradiated at 50°C inside a uranium tube in the 
NRX reactor (Cook and Cushing 1953). Half the fast neutrons hitting the 
sample were from the uranium tube and had a fission energy spectrum. 
The other half from neighbouring fuel rods had the ususal E~ spectrum. 
The definition of fast neutrons used in this report is neutrons with energy 
above 1200ev. This is the minimum energy neutron capable of knocking 
a platinum atom from its lattice site. 

After irradiation, the platinum wires were decontaminated and platinum 
potential leads spot-welded on. 

02 
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§ 4. Mernops ror Dretectine DEFECTS IN PLATINUM 
4.1. Electrical Resistivity 


Potential leads were spot-welded to the wire samples, which were 
0-010 in. in diameter, and the electrical resistances measured in liquid 
nitrogen. The values obtained were corrected to a constant nitrogen bath 
temperature by measuring the resistance of an annealed sample in the 
same bath. The electrical circuit was made completely of copper with 
low thermal e.m.f. solder at all joints. A 100ma current was used with 
the bath vigorously stirred to prevent temperature fluctuations due to 
gas bubbles forming on the sample. The potential was read to 0-1 
microvolt. By reversing the sample current and potentiometer battery 
to eliminate stray e.m.f.’s, the resistance values were reproduceable to 
one part in thirty thousand. 

The samples were annealed in silicone oil up to 250°c and in molten salt 
between 250° and 900°c. 


‘NORMAL 
BALANCE 
ENCLOSURE 


PLASTIC 
ENCLOSURE 


LIQUID LEVEL 


Schematic diagram showing the balance assembly for density measurements. 


4.2. X-ray Line Width 
A long 0-010 in. diameter platinum wire was extended 8°% in liquid 
nitrogen. Sections of it were spot-welded side by side on a platinum 
frame to fit the diffractometer. The half widths of the (331), (420) and 
(422) x-ray lines were then measured after annealing the wire sample 
assembly at successively higher temperatures. 


4.3. Density 
The density of platinum was determined by weighing a 13g sample in 
air and water with a weight sensitivity of 0-01 milligrams. The density 
of the sample was directly compared with an annealed reference sample of 
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the same size using the system shown in fig. 3. The sample was put in the 
upper left pan and the reference sample in the upper right pan to obtain 
the difference between their weights in air. The sample was then placed 
under the liquid in the lower left pan and the standard in the lower right 
pan, to obtain the difference between their weights in liquid. The cal- 
culations for this method of density determination are given elsewhere 
(Piercy 1958). This is a very sensitive method for measuring changes in 
density since the measurements are not strongly dependent on the liquid 
density, air density, or temperature. Of course, the absolute value 
depends on the accuracy of the density taken for the reference sample. 

Distilled water with the air removed was used for the density measure- 
ments. However, hysteresis in the contact angle between water and the 
0-003in. tantalum wire, caused the measured weight to depend on the 
direction of vertical motion of the balance pan. This effect was markedly 
reduced by using a solution of 0-08% wetting agent (Kodak Photoflo) in 
distilled water, which resulted in a density reproducibility of 0-005% for a 
13g platinum sample. The addition of the wetting agent did not affect 
the absolute value within this limit of error. 


4.4. Reaction Kinetics 

Once the excess point defects were introduced into the sample, the 
preliminary investigation of measuring the drop in resistivity due to 
annealing the sample for equal times at successively higher temperatures 
(an isochronal anneal) determined the minimum temperatures at which 
these defects were mobile. The decrease in resistance for each 5min 
anneal was plotted as a function of the annealing temperature, a typical 
result for deformed platinum being shown in fig. 4. Once the tempera- 
ture range for recovery was known, the activation energy was found by 
measuring the change in slope of the isothermal resistivity—time curve 
when the temperature was raised by 10°c. The value for the activation 
energy AH, was calculated from the usual equation 


Lek 
A= Fr, 7, 108. (82/1) . . . . . (1) 


_where K is Boltzmann’s constant, S, is the final slope of the resistivity— 
time curve at absolute temperature 7', and S, is the initial slope after 
raising the temperature from 7’, to T,. This process was repeated until 
there was no further drop in resistivity. A typical set of curves for 
irradiated platinum is shown in fig. 5. 

The slope at the end of each curve was determined in the following 
manner: the slope at the midpoint of a curve between two experimental 
points A and B was assumed equal to that of the straight line joining A and 
B. Therefore by measuring seven experimental points on. the isothermal 
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Isochronal recovery of platinum after deformation in liquid nitrogen. The 
drop in sample resistance for successive anneals of 5 min at temperature 
T is plotted against annealing temperature. 


Fig. 5 
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A typical set of curves to determine the activation energy required for recovery 
of irradiated platinum. 
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anneal, the values for the slope of the curve at the six intermediate 
points were determined. Then, by plotting the curve for these values of 
slope against the resistivity at these intermediate points and extrapolating 
the curve a short distance, the value for the slope at the end of the curve 
was obtained. The logarithm of the slope against the resistance at the 
intermediate points was usually plotted in order to have a straighter line 
to extrapolate. A typical example for irradiated platinum is shown in 
fig. 6. This is one of the few methods that determines the slope using all 
the experimental points equally. 


SR/OT MICRO-OHMS PER SECOND 


0.0206 0.0208 0.0210 
RESISTANCE IN OHMS 


Determination of slope at the end of an isothermal annealing curve. The log 
of the slope found at a point midway between two experimental points 
is plotted against the electrical resistance at that point. The line is 
extrapolated to point C to obtain the slope at the end of the isothermal 
curve. 


When the activation energy was known, the number of jumps required 
by the defect to reach the sink was determined from the slope of the curve 
obtained by plotting the logarithm of the fraction of resistivity increase 
remaining against the time at constant annealing | temperature and 
applying eqn. (9) of the Appendix. The order of reaction was found from 
the same isothermal anneal using the method of Meechan and Brinkman 


(1956). 
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§ 5. RESULTS AND OBSERVATIONS 
5.1. Changes in Physical Properties 
The electrical resistivity isochronal recovery curve for platinum ex- 
tended 8% in liquid nitrogen (fig. 4) has two distinct recovery stages; one 
below 100°c in which 0-004 micro-ohm-cm recovers, and one between 
100°c and 600°c in which 0-024 micro-ohm-cm recovers. For comparison, 
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Recovery of electrical resistivity, hardness and x-ray line width of platinum 
samples deformed 8% in liquid nitrogen. Successive anneals for 30 min 
every 100°C were used. 


the electrical resistivity of platinum in liquid nitrogen is 1-85 micro-ohm- 
cm. There was no appreciable recovery below —50°c which is contrary 
to the experiments of Manintveld (1954). The recovery of x-ray line 
width, electrical resistivity and hardness for platinum deformed 8% in 
liquid nitrogen are shown together in fig. 7. Because the temperature 
intervals are larger, this curve for the recovery of electrical resistivity does 
not resolve the two peaks shown in fig. 4. The x-ray line width does not 
change until nearly all the increase in electrical resistivity has recovered 
indicating that both stages in the recovery of electrical resistivity shown 
in fig. 4 are mainly due to point defects, rather than to annihilation or 
redistribution of dislocations. However, some annihilation of dislocations 
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may be occurring at the very end of the recovery of electrical resistivity. 
The recovery of the hardness and x-ray line width above 400°c is probably 
due to recrystallization. 
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Isochronal recovery of quenched platinum showing the drop in sample resistance 
for successive anneals of 15 sec at temperature 7’. 


The isochronal recovery of quenched platinum shown, in fig. 8 has one 
recovery stage in the region 300 to 600°c, during which a resistivity drop 
of 0:12 micro-ohm-cm occurs. There was no further recovery on annealing 
up to 800°c. The isochronal anneal of irradiated platinum shown in 
fig. 9 also has only one recovery state and this occurs in the region 130- 
350°c. However, the large decrease of 0-48 micro-ohm-cm occurs here. 
There was no further recovery on annealing up to 900°c. 
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The recovery of the density and yield stress of irradiated platinum was 
examined to determine whether these recovered at the same rate as the 
electrical resistance. The results shown in table 1 indicate that both the 
density and the yield stress recover at the same rate as the electrical 
resistance and, therefore, are probably due to the same defect. Since the 


Fig. 9 


{000 


‘ 


@ 
ce) 
Oo 


8 


400 


200 


RESISTANCE DROP IN MICRO-OHMS DURING IO MINUTE ANNEAL AT TEMPERATURE T 


100 200 300 
TEMPERATURE T IN °C 


Isochronal recovery of irradiated platinum showing the drop in sample resistance 
for successive anneals of 10 min at temperature 7’. 


density decrease is 0-07°% due to irradiation, while the lattice parameter 
change for the same irradiation is less than 0:005% (Tuxworth 1958 
unpublished data), the defects remaining in the sample are probably single 
double, or clusters of vacancies. 
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From the results shown in table 1, the ratio of the resistivity change to 
fractional volume change due to irradiation is 9x 10-*ohm-cm. Since 
this is comparable with the value 3x 10~*ohm-cm found in gold wire 
that was quenched from 900°c (Bauerle and Koehler 1957), the defects in 
irradiated platinum are probably similar to those in quenched gold, i.e. 
isolated or small grouping of vacancies rather than large clusters. If one 
vacancy in the irradiated platinum produces an effective volume change 
of 0-5 atomic volumes (the number will be less than 1-0 due to relaxation 
about the vacancy), then the resistivity per atomic per cent vacancies 
when calculated from the ratio of resistivity change to fractional volume 
change is 4-5micro-ohm-cm. This value, which is about three times that 
predicted for copper (Overhauser and Gorman 1956) is a reasonable 
resistivity for vacancies in platinum and will be used for later calculations 
of defect concentrations produced in platinum by other means. The 
concentration of vacancies determined for the irradiated sample by using 
this value is 1-2 x 10-3. 

The defect concentration produced by fast neutrons can be calculated 
from the theory of Seitz and Koehler (1956), and also the theory of 
Brinkman (1956). Taking the total elastic scattering cross section for 
fast neutrons to be 5:5 barns (Langsdorf et al. 1957), the concentration of 
defects produced in platinum which is irradiated with 2~x 107° fast 
neutrons per square centimetre in NRX is 0-30 when calculated from the 
theory of Seitz and Koehler. On the other hand, Brinkman’s estimate — 
for the concentration, of defects remaining in a displacement spike is only 
about 0-001 (Brinkman 1954). Therefore, this should be the limiting value 
when all the atoms of the sample have been in at least one displacement 
spike, as would probably occur with the high integrated flux used here. 
The concentration of defects measured in this experiment, however, is 
expected to be much lower than the actual numbers produced, due to 
annihilation of the vacancies by interstitial atoms that are mobile during 
irradiation. Therefore, the experimental value of 0-0012 deduced from 
the density and resistivity change indicates that Brinkman’s estimate of 
0-001 for the concentration of defects remaining in a displacement spike 
is much too low if the displacement spike concept is applicable to platinum. 


5.2. Reaction Kinetics 


The results for the activation energy measurements on extended, 
quenched and irradiated platinum are given in table 2. They were all 
measured by the change in slope of an isothermal recovery curve when the 
temperature was raised. 

A lattice distortion or strain will alter the distance between atoms and 
hence may alter the activation energy required for the atomic jump 
process. Since a comparison between the activation energy of recovery 
after cold working, quenching, and irradiating is required, and these do 
not produce the same lattice strain in the metal, the possible change in 
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activation energy with lattice strain was investigated using an irradiated 
sample. A load equal to half the proportional limit of platinum at 400°c 
was hung from an irradiated sample during recovery. Although this 
increased the rate of recovery by a factor of two, the activation energy was 
still the same within experimental error. 


Table 2. Activation Energy for the Recovery of Platinum after 
Deformation, Quenching or Irradiation 


TPecectricnt Temperature Activation Number of 
range °C energy ev | determinations 
Extended 8—10°% in liquid 
nitrogen 38 to 135 0-73 + 0-067 5 
Extended 8-10°% in liquid 
nitrogen and annealed to 
remove lower recovery 
stage 230 to 360 1-43 + 0-09 7 
Quenched from above 
1600°c to 20°c in less 
than 0:04 sec 240 to 370 1-13+0-11 14 
Trradiated at 50°c with 
6x10!% fast neutrons 
per square centimetre 190 to 280 1-40+0-11 bt 


+ The limits given are the standard deviations of all measurements taken. 


The average values of activation energy from table 2 were used to 
calculate the order of reaction and the number of jumps required by the 
defect to reach the sink. These, together with the concentration of defects 
‘(based on recovery of electrical resistivity, assuming 4:5 micro-ohm-cm 
per atomic percent defects) are shown in table 3. The activation energy 
for the recovery of irradiated platinum, the only case for which two 
identical samples were available, was also measured by the technique of 
Meechan and Brinkman (1956). The result shown in fig. 10 and given in 
table 3 is within the limit of error of the slope change method. 


§ 6. DIscUSSION 


A comparison between the present results for the reaction kinetics and 
other investigators are given in table 4. The values for activation energy 
do not agree except for quenched platinum. Since the purity of 
Manintveld’s extended platinum was not reported, the difference in results 
may be due to an impurity effect. Unfortunately, his annealing experi- 
ments were not taken above 90°c. Dugdale’s results on platinum irradi- 
ated or extended at room temperature (Dugdale 1952) are both the same 
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Fig. 10 
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Activation energy for the recovery of irradiated platinum determined by the 
method of Meechan and Brinkman (Meechan 1956). The time Ar 
required at 196°c for the same recovery occurring 10 min at temperature 
T of the isochronal anneal is plotted against 1/7’. The slope of the line 
gives an activation energy of 1-46 ev for the whole recovery process. 


but are significantly below those reported here. However, he made only 
three activation energy measurements and all were taken after about 8% 
recovery. Once again, the difference may be due to sample impurity since 
his is unspecified. Because of these discrepancies, no single model can 
explain all the past results. Therefore a model will be proposed to explain 
the results of the present experiments in which point defects were produced 
by three methods and in which all measurements were made under the 


same conditions. 
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Since the density of platinum decreased after irradiation while the lattice 
parameter showed no significant change, the defects remaining after a 
room temperature irradiation are probably vacancies, divacancies or 
vacancy clusters. Because the changes in density and _ electrical 
resistance recover at the same rate after irradiation, they are probably 
caused by the same point defect. Therefore, the recovery of electrical 
resistivity in irradiated platinum that requires an activation energy of 
1-43 ev is probably due to vacancies, divacancies, or vacancy clusters. 

There are several arguments that strongly suggest the defects in irradi- 
ated platinum are isolated or small groups of vacancies rather than large 
clusters. The first, mentioned previously, is the agreement between the 


Table 4. A Summary of Activation Energies Measured for the Recovery 
of Platinum 


Recovery 
Method of temperature | Activation 
introducing range energy Reference 
defects examined ev 
¢ 
Extended at —196°c — 150 to —70 0-22 Manintveld 1954 
—10 to 70 0-99 Bs 
Extended at 20°c 70 to 100 1-19 Dugdale 1952 
Extended at —196°c 0 to 100 0-73 Present work 
200 to 400 1-43 a Bs 
Quenching from above — LI Lazarev and 
1600°c Ovcharenko 1955 
300 to 500 1-10 Bradshaw and Pearson 
1956 
300 to 500 1-13 Present work 
Quenching from 1100°c 450 to 580 1-42 Ascoli et al. 1958 
to 1400°c — 1-48 Bacchella et al. 1959 
Reactor irradiation at 70 to 90 In19 Dugdale 1952 
50°C 100 to 300 1-43 Present work 


ratio of density decrease to resistivity increase in irradiated platinum 
with that found in quenched gold. The second, is the excellent agreement 
between the migration energy for defects produced by irradiation (1-43 ev) 
with those produced by quenching platinum from temperatures substanti- 
ally below the melting point, 1-42 ev (Ascoli et al. 1958). A quench from 
these temperatures is expected to leave single vacancies while a quench 
from near the melting point may leave double vacancies (Koehler et al. 
1957). The lower values of 1-lev required for defect migration after 
quenching from near the melting point as shown in table 4 are probably 
due to double vacancies. ; 
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The third argument that suggests the defects produced by irradiation 
are single vacancies, is the qualitative agreement between the measured 
migration energy of 1:43+0-lev and the calculated migration energy of 
1-65 + 0-15 ev from the energy of self-diffusion in platinum of 2-96 + 0-06 ev 
(Kidson and Ross 1957) less the formation energy of vacancies in platinum 
1-30 +0-lev (Bradshaw and Pearson 1956, Ascoli et al. 1958, Bacchella 
et al. 1959). 

The activation energy required for recovery of the electrical resistivity 
after 8% extension in liquid nitrogen was 0-73 and 1:43ev. The higher 
value has already been attributed to the migration of single vacancies and 
the lower one is perhaps due to the release of interstitial atoms from trap- 
ping sites. Both vacancies and interstitial atoms may be produced by the 
crossing of screw dislocations during deformation. 

The order of reaction is a measure of the rate of change of the number of 
jumps required by a defect to reach the sink during the recovery process. 
Once it rises above the value two, its interpretation becomes vague and 
in many cases may be caused by a nonrandom initial distribution of defects. 
This nonrandom distribution of defects is one possible explanation, for the 
high order of reaction and the low number of jumps required for the 
recovery of irradiated platinum. If the damage from a fast neutron is 
concentrated in one small volume the local concentration of defects may be 
very high. Therefore the vacancies would require only a small number of 
jumps before combining with another vacancy, after which they could 
move more rapidly. Also, if the damage from fast neutrons is concentra- 
ted in a small volume, there may be other traps for vacancies such as 
dislocation rings or clusters of interstitial atoms that are produced by 
the same neutron. Therefore the vacancies would again require only a 
small number of jumps before reaching these sinks. 

From the number of jumps required for the defect to reach a sink in the 
recovery of quenched and deformed platinum, an interpretation of the 
type of sink can be made. In the quenched sample with a defect concen- 
tration of 3 x 10-4, the number of jumps (8 x 108) is much too high for 
recombination of defects. The low order of reaction 1:46 also dis- 
courages this interpretation. However, if the defects are going to dis- 
locations, an argument which agrees with the observed number of jumps, 
an order of reaction equal to 1-0 might be expected. The higher order of 
reaction observed (1:46) may be caused by the change in dislocation net- 
work due to extensive climb of dislocations which absorb defects. Thus, 
the last defects will migrate to a much different distribution and perhaps 
even different type of dislocation network than was initially present. 

The eight per cent extension of platinum is expected to increase the 
dislocation density from 108 to about 101° lines per square centimetre. 
Therefore, the number of jumps required by the defect to reach a dis- 
location will be reduced from 107 to 10° (Lomer and Cottrell 1955) which is 
similar to the number of jumps required for the defects to combine. 
The recovery of extended platinum at 60°c appears to be by defects 
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combining since the order of reaction is close to two. However, the high 
order of reaction (5-05) for the recovery of deformed platinum at 300°c, 
which is reflected in the large increase in number of jumps during recovery, 
indicates the defects are diffusing to several different types of exhaustible 
sinks or the sink concentration is decreasing during the recovery. The 
latter view is perhaps more favourable, since the experiment on the 
recovery of x-ray line broadening shown in fig. 7 indicated that some 
reduction in, the dislocation density may occur during the recovery of the 
electrical resistivity. However, because of the unique activation energy 
and the large drop in line broadening on annealing at higher temperatures, 
the main recovery of electrical resistivity in the region of 300°c is thought 
to be due to vacancies rather than to dislocations. 


§ 7. CONCLUSIONS 


An interpretation of the point defects in platinum consistent with all 
our experimental results is the following : 


(1) The principal lattice defects remaining in platinum after fast 
neutron irradiation are single vacancies with an activation energy for 
migration of 1-43 ev. 

(2) The principal lattice defects remaining in platinum after quenching 
from above 1600°c are divacancies with an activation energy for migration 
of 1-13ev. 

(3) The two principal lattice defects remaining in platinum after an 
extension of 8% in liquid nitrogen are an unknown defect, perhaps inter- 
stitial atoms, with an activation energy for migration of 0-:73ev, and 
single vacancies with an activation energy for migration of 1-43 ev. 
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APPENDIX 
NuMBER OF JumpPsS REQUIRED BY DEFECT TO REacH SINK 


The average number of jumps required by defects to reach a sink can 
be calculated from the recovery of electrical resistivity at constant tem- 
perature 7’. If the defect concentration is assumed proportional to the 
additional resistivity R above that for the fully recovered sample, the 
recovery rate depends on three factors 


dRidi= EP pe Se) J ee 
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where ¢ is the time, I’ the defect jump frequency, and p the probability 
per jump that a defect reaches a sink. 


The jump frequency IT also depends on three factors 
die ANY) OX ( NM lyme eas (DQ) 


where Z is the number of adjacent sites to which the defect may go, 
av is the frequency of attempts by the defect to move to each adjacent 
site, and exp (— AF,,/K 7’) is the fraction of time that a defect has sufficient 
free energy to jump. =v is the vibration frequency of the atom, and « is 
the fraction of these vibrations in each jump direction. The free energy 
required by the defect to jump is 


ARR TAS We ee a a ees) 


where AH, is the activation energy required to move the defect one jump 
and AS, is the entropy change produced by the defect moving from one 
site to the saddle point midway between sites. 

The probability per jump that a defect reaches a sink is equal to the 
fraction, of lattice sites from which the defect is directly attracted to the 
sink. This depends on the concentration of sink sites, their distribution, 
and the maximum distance of attraction between the defect and sink. 
Therefore 


P=) *=BCs ee eet wt hte tee) 
where j is the average number of jumps required by the defect to reach the 
sink, C's is the concentration of sink sites and f is a term depending on the 
sink distribution and the maximum distance of attraction between defect 
and sink. The term BC, could well be called the ‘effective sink concen- 
tration ’. 

Putting eqns. (2), (3), (4) into eqn. (1) gives: 
g 1dk__d(mk) 


=Zavj—texp (AS,/K)exp(—AH,/KT) (5) 


R dt dt 
where S is the slope of the graph In Ff vs. time. Solving for 7 gives: 
4 =Zav S exp (AS,/K) exp (— AH,/ KT). Ra ac a.) 


AH, is calculated from other experiments using eqn. (1) of the main text. 
Unfortunately the frequency factor Zav exp (AS,/K) has not been measured 
directly, but each factor can be estimated. Because the total vibrational — 
energy is considered in the AH, term in eqn. (6), v must be the frequency 
of the resultant vibration containing the total vibrational energy of the 
atom making the jump and « is the fraction of these resultant vibrations in 
each jump direction. The Debye frequency is usually used for v. Bradshaw 
and Pearson (1956) and LeClaire (1953) take «=1 but the value a=Z™ 
assumed by Jost (1952) seems more reasonable for vacancy diffusion since 
only the resultant vibration is considered. Because there are Z possible 
jump directions, the fraction of the resultant vibrations in any one jump 
direction will be about Z-!. If the resultant vibration of an atom adjacent 
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to a vacancy is always in the direction of the vacancy, then «=I. This is 
unlikely, however, because each atom vibrates in phase with its nearest 
neighbours and the wavelength of vibration is usually greater than three 
to four atoms. The remaining factor, exp(AS,/K), was assumed by 
Bradshaw and Pearson (1956) to lie between one and ten. 

Because of uncertainty in the choice of these parameters for the fre- 
quency factor, a more direct value was taken from the self diffusion 
measurements of Kidson and Ross (1957). Their measured value for Dy 
was 0:-33cm2sec"! which gives (LeClaire 1953) 


avexp [(AS,+ AS,)/K]=D,/a? = 2-2 x 104 es co ade came tad 


where a= 3-9 x 10-8 cm is the lattice parameter and AS, is the vibrational 
entropy. due to a vacancy at a lattice site. The left side of eqn. (7) is 
similar to the frequency factor in eqn. (6) except for the added AS, term. 
However, from calculations on the relative values of AS, and AS, 
(Huntington et al. 1955) 


exp (AS,/K) =0-24exp [(AS, + AS,)/K]. + tala te SD 


Therefore by combining eqns. (7) and (8) with (6) and taking Z= 12 for 
f.c.c. lattices, 7 becomes 


j=6-3 x 104 S-texp(<AH/KT) . .°. . . (9) 


where S is the slope of In # vs. time and AH, is the activation energy for 
recovery at temperature 7’. This equation was used to calculate the 
number of jumps required by the defect to reach the sink for all recovery 
stages investigated. The frequency factor may be incorrect if recovery 
other than vacancy diffusion occurs. This frequency factor is the same 
as that of Bradshaw and Pearson’s (1956) using their assumptions 
(B=exp (AS,/K)=5, «=1, v=10'%) and including the extra Z term 
discussed in his later paper (Bradshaw and Pearson 1957). 

An estimate of B=exp (AS,/K) can be made from eqns. (7) and (8). 
If «=Z, which seems reasonable for vacancy diffusion, and v=5 x 1022, 
the Debye frequency for platinum (Seitz 1940); then a value for B of 125 
is required to match the diffusion results. 
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ABSTRACT 


Diffraction patterns of low energy electrons produced by a (111) face of 
a crystal of nickel have been accelerated sufficiently to show on a fluorescent 
screen, where they can be photographed. Various gases at low controlled 
pressures have been admitted to a clean surface, and the arrangements of 
atoms or molecules upon the surface deduced from their diffraction patterns. 
Molecules of CO are adsorbed at low pressures as monolayers. They form 
a single two-dimensional crystal having a spacing between molecules double 
that between surface nickel atoms ; these molecules are completely removed 
from the surface by heating to 175°c.. Experiments with nitrogen have 
produced a polycrystalline monolayer having the same structure as the CO 
crystal. (Because of low sticking probability adequate experiments to prove 
that the layer is really nitrogen are difficult and have not yet been 
satisfactorily completed. A new type of experimental equipment seems 
to be required which is now being constructed.) Crystal growth in this 
monolayer is induced by heating to 125°c; heating above 200°c results 
in transformation into one or the other of two polymorphic forms either 
of which is more stable than the original structure ; heating to 300°c 
results in melting the monolayer without removing it, the crystalline 
structure returning upon cooling. Oxygen can be adsorbed in a layer 
several atoms or molecules thick ; this shows no structure until annealed 
at 125°c after which the structure is that of the most complex of the 
polymorphic forms assigned to nitrogen; this complex structure is 
polycrystalline with very small individual crystals ; its diffraction pattern 
is unchanged when the oxygen is reduced to a single monolayer by heating. 
No diffraction pattern has been obtained from hydrogen, although the 
presence of hydrogen is proved by rise in pressure when the crystal is 
subsequently heated. 


§ 1. LyrropuctTIon 


THE arrangement of foreign atoms adsorbed as a single monolayer upon 
a crystal surface can be determined by the diffraction of very low energy 
This has been known for a long time 
(Davisson and Germer 1927, Germer 1929) and a great deal of work in this 
field has been carried out by Farnsworth et al. (1929 to 1959). 
present time all observations of this sort have been obtained by the 
laborious method of discovering the features of a diffraction pattern by 
capturing scattered electrons in a Faraday box which is moved about over 


the entire field of the pattern. 


+ Communicated by the Authors. 
} Now at Georgia Institute of Technology, Atlanta, Georgia. 


Up to the 
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Ehrenberg (1934) showed that diffraction patterns produced by low 
energy electrons can be accelerated through a considerable potential 
difference so that the entire pattern can be seen on a fluorescent screen. 
This represented potentially a very great improvement in the experimental 
procedure over the earlier method. Nevertheless, Ehrenberg did not 
develop the technique to a point where it had practical usefulness, and he 
did not discover patterns due to adsorbed atoms. 

We have improved this means of observing a diffraction pattern of 
low-energy electrons sufficiently to make it a valuable research tool, and 
by means of it we have determined the arrangements of various kinds of 
atoms and molecules upon a crystal surface. With this new technique, 
experiments which were very difficult by the older method are now quite 
easy, and a vast number of experiments can be carried out which were 
impossible up to this time (American Crystallographic Association 1959). 


Fig. 1 
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Sketch of the arrangement of the electron diffraction tube, drawn to scale. 


Details of the experimental tube will be published elsewhere. or the 
present purpose, the sketch of fig. 1 will suffice. Electrons emerging from 
the gun strike the crystal at normal incidence. The diffraction pattern 
moves in a field free space until it reaches a fine-mesh grid. Between this 
grid and a second grid is a decelerating field which removes most of the 
slow electrons. Beyond the second grid the pattern is accelerated by 1000 
or 2000 v to a fluorescent screen. The pattern can be observed directly on 
the screen, or by its reflection in the mirror which is inclined by 45° to the 
direction of the primary beam. The reflection can be photographed. 

The experimental tube of fig. 1 is attached to vacuum pumps and to a 
Granville—-Phillips type vacuum valve through which various gases can 
be admitted at low controlled pressures. The crystal can be outgassed 
by electron bombardment and the pressure of residual gas reduced to 
about 2x 101mm Hg. Diffraction patterns can be observed when the 
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surface of the crystal is completely clean and the gradual development of 
new patterns noted as atoms or molecules of any one of a number of gases 
settle upon the surface The observation of these changes is made easy 
by the fact that the separations of adsorbed gas atoms is greater than the 
separations of nickel atoms in the surface layer, and, as a consequence of 
this, diffraction patterns from adsorbed gas atoms are observed in general 
in a range of voltages lower than the lowest voltages for which the nickel 
lattice can give any pattern. Thus in most cases there is no superposition 
of a diffraction pattern from adsorbed gas atoms upon a pattern arising 
from the nickel crystal and hence no possible confusion. 


§ 2. DirFRACTION PATTERNS FROM A CLEAN CRYSTAL 


Representative diffraction patterns from a clean (111) face of a nickel 
crystatare reproduced in fig. 2 (Pl. 26). The first three of these, fig. 2(A), 
(B) and (C), represent respectively (533), (644) and (466) reflections from 
the lattice of the crystal, the patterns being taken at the electron beam 
voltages, 129, 200, and 262, for which the diffraction spots are strongest. 
Each pattern has threefold symmetry (except for some skewness due to 
an imperfect electron beam, and the fact that the spot at about 10 o’clock 
in fig. 2(C) is obscured by the position of the electron gun). In each case 
the diffraction spots appear in one or the other of the two principal azimuths 
of the crystal, designated A and B in the sketch of fig. 3. These are the 
azimuth designations used in earlier experiments (Davisson 1928a,b, 
1929, Germer 1928). The A azimuth has also been called the (111) azimuth, 
and the B azimuth the (100) azimuth (Davisson and Germer 1927). 

In each of these three patterns, weak diffraction spots are seen in the 
other principal azimuth. This occurs because of the poor resolving power 
due to limited penetration of the electron beam into the crystal. This 
effect is shown in exaggerated form in fig. 2(D) in which the diffraction 
spots appear with equal intensities in the two azimuths; the beam voltage 
in this case is intermediate between 106 v, at which (244) reflections appear 
at maximum intensities in the A azimuth, and 129v, at which (533) 
reflections appear at maximum intensities in the B azimuth. 

For each diffraction pattern the tangent of the angle ¢ between the 
primary beam and a diffraction beam is obtained by dividing the distance 
of a diffraction spot from the center of the pattern (measured on such a 
photograph as those of fig. 2) by 44mm, the distance from the crystal 
face to the fine-mesh grid, and multiplying by the photographic enlarge- 
ment of the pattern. The latter is the ratio of the diameter of the grid 
measured on the photograph (easily seen on fig. 2(D)) to 67-5 mm, the true 
diameter of the grid. 

The lowest voltage diffraction beams from a clean crystal appear at 
maximum intensity at 106v, in the A azimuth. When the voltage of the 
primary beam, JV, is increased steadily from this value to 400, the limit of 
our observations, sets of three diffraction beams appear alternatively in 
the two principal azimuths, each set increasing gradually in intensity to a 
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Plots of all diffraction beams observed from the nickel crystal. 
iffraction beams plotted as crosses ; ca 
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maximum value then decreasing in intensity and disappearing. For each 
such set of beams the electron wave-length at maximum intensity, cal- 
culated from A= (150/V)?, has been plotted in fig. 4 against the sine of the 
co-latitude angle, ¢, measured as described above. The observed dif- 
fraction beams at maximum intensity are plotted as crosses and the 
calculated positions of all possible Laue x-ray diffraction beams in this 
wavelength range are plotted as open circles. The region of the plot 
which is inaccessible to observation is marked off by cross hatching. The 
lines in fig. 4 upon, which the experimental points fall have a slope of 2-154 
which is the spacing of lines of surface atoms in the A and B azimuths. 

The correct correlation between observed and calculated diffraction 
beams is quite clear. All possible diffraction beams are found. These 
are written down in table 1, with the potentials at which they are observed 
at maximum intensity, and the theoretical potentials for zero value of the 
inner work function. The experimental values of inner potential probably 
do not differ from the average value of 16v by more than the possible 
experimental errors of finding the potential for maximum intensity. Thus 
the experiments supply no data on variation of inner potential with angle 
or voltage. 


Table 1. Laue Diffraction Beams from a Nickel Crystal (See fig. 4) 


Azimuth A Azimuth B. 
Voltage Voltage 
“tion | 0% Cale. | potemtiat | tion | 0% Cale. | potential 
244 106 ALo 13 533 129 140 1] 
355 178 188 10 644 200 215 15 
466 262 276 14 755 294 310 16 
577 355 = 3383 28 866 403 424 21 


Average inner potential 16 v. 


§ 3. ADSORBED LAYERS 


The nickel crystal had received considerable heating before the diffraction 
patterns of fig. 2 were obtained. A very brief summary of this heating 
seems to be a necessary prelude to an account of the diffraction patterns 
obtained from gas films deposited upon a clean surface. 

Before the crystal was heated at all, no diffraction pattern could be 
observed. After heating to only 500°C, all of the Laue diffraction beams 
listed above were found, but there were in addition very elaborate patterns 
at low voltages. An example of one of these patterns is reproduced as 
fig. 5 (Pl. 28). (The structure giving rise to these complex patterns was 
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later identified as having the unit cell of the ‘ Q-structure’ which is shown 
in fig. 6). Low voltage diffraction patterns were entirely eliminated only 
after the crystal had been heated a number of times at maximum tem- 
peratures of 830°c. There were finally no traces of patterns other than 
the normal Laue beams from the nickel crystal, not even a ‘muddy’ 
appearance of an annular region of the fluorescent screen, appearing at 
low voltages, which was later associated with less than a monolayer of 
gas on the surface. 


Fig. 6 


4 3 


Sketches of the three different arrangements which atoms, which are presumably 
nitrogen, can take up upon the (111) surface of a nickel crystal. Black 
dots represent surface nickel atoms, open circles adsorbed atoms. 
The unit cells of the three different structures are drawn in; in the 
case of the Q structure, locations of atoms within the unit cel) are not 
known. 


3.1. Nitrogen on Nickel 


Nitrogen was the first gas which was admitted to the experimental tube 
with the nickel crystal in a condition believed to be clean. Three different 
types of low voltage diffraction patterns have been observed after the 
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admission of nitrogen, corresponding to three different arrangements of 
atoms in a monolayer on the surface. One of these structures is not yet 
completely understood, although its two-dimensional unit cell is known. 
Sketches of these structures are shown in fig. 6. 

Because of the rather low sticking probability observed in these experi- 
ments it has not been established with certainty that the adsorbed gas is 
nitrogen rather than an impurity gas present in small amount but having 
sticking probability close to unity. This uncertainty is aggravated by 
uncontrolled fluctuations of pressure incidental to. the method of cleaning 
the crystal at high temperatures. The bombarding tungsten ribbon 
adsorbs gas readily, and the release of this gas whenever the ribbon is 
heated raises the pressure to undesirably high values. A nickel crystal 
which can be heated by conduction will reduce this difficulty, and this seems 
to be essential for better control of gas pressure. 

In all of the following experiments with nitrogen the adsorbed gas will 
be spoken of as nitrogen, which we now doubt. Presumptive evidence that 
it is atomic rather than molecular is furnished by the low sticking probability 
and by the fact that it is not removed until crystal temperatures above 
400°c are reached. 

In the sketches of fig. 6 the black dots represent the centers of nickel 
atoms lying in the surface layer, and the open circles the positions of 
adsorbed nitrogen atoms. In the structure marked 4 the separations 
between adjacent nitrogen atoms is double that between, nickel atoms in 
the surface layer ; the ratio of the number of nickel atoms to the number 
of nitrogen atoms is 4, which number has been chosen to designate the 
structure. In the structure marked 3 the corresponding ratio of distances 
is 342, and the ratio of the numbers of atoms is 3. 

The structure marked Q is much more complex with a two-dimensional 
hexagonal unit cell having an edge of 12-44, five times larger than the 
separation of adjacent nickel atoms. The relative intensities of different 
diffraction spots of a pattern representing the Q structure are not the same 
in different experiments. For example, the (22) and (33) reflections are 
sometimes strong and sometimes weak. This seems to mean that the 
locations of atoms within the unit cell of this structure is not always the 
same. Furthermore, there is, in general, no single arrangement in, the 
different unit cells which can account for the observed intensities, which is 
interpreted to mean that the arrangement is not the same in all the unit 
cells. 

Tn, all three of these structures it has been assumed that adsorbed atoms 
are located between surface nickel atoms as indicated, although the 
diffraction patterns do not establish this fact. As the sketches of fig. 6 
are reproduced here, there is no distinction between the centres of small 
triangles of nickel atoms having their apices up and down. If we identify 
azimuths in fig. 6 with those marked on fig. 3 (and with the approximate 
orientation of the diffraction patterns of fig. 2(A), (B) and (C), then, the 
placing of atoms at the centres of triangles with apices upward (as is done) 
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is equivalent to placing them above nickel atoms of the second layer below 
the surface rather than above atoms only one layer below. There is some 
experimental evidence that this is the correct placement. 

Representative diffraction patterns from which the 4, 3 and Q structures 
have been deduced are reproduced respectively in figs. 7, 8 and 9 
(Pls. 27-29). Data relating to the patterns representing the 4 and 3 
structures are assembled in table 2, and data relating to the patterns 
from the Q structure in table 3. 

The experimental conditions under which atoms take up these different 
arrangements on the crystal surface are imperfectly understood, and the 
clarification must wait upon further work. The observations which are 
now available are summarized below, together with the difficulties 
encountered in understanding the factors involved. 


Table 2. Patterns Representing the 4 and 3 Structures 


(1) (2) (3) (4) (5) (6) (7) (8) 


Voltage pane Seer aoe Observed beams 
hanes Observ - 
length eee eres Spacing ~Azimuth}| Spacing Azimuth 
The 4—structure 
5o—y From ,2 
Fig. 7(A) 1-69 4 6-84 (10) 4-304 ASB 4-304 ASB 
. to 2-94 
ae From 
. 57 Vv 5 : an Wie esiegolaves 
Fig. 7(B) 1-615 4 eee (10) 4-304 A,B 44580 Fe cay 
0 7: 
From : 
=. 65 Vv : Ann “Gm i B strong 
Fig. 7(C) L524 on, (10) 4-30 A A, B ae Oe Nears 
0 2-6 
77 From ; 
Fig. 7(D) Lee 5-7 A (10) 4:304 A,B 4-514 A,B 
oe 2-44 11 2-484 C 2-464 
1-404 dk (11) C 
The 3—structure 
62 v From A ; ne ’ 
ij re “BA 1¢ 3-724 y 3-764 } 
Fig. 8 15554 6-34 (10) 7 


The 4 structure seems to be favoured by deposition at very low pressure 
upon a crystal at room temperature. At an appropriate electron beam 
voltage, the pattern representing this structure makes its appearance on. 
the fluorescent screen rather gradually. The diffraction spots are initially 
very diffuse, presumably representing simultaneous nucleation of individual 
two-dimensional crystallites at many different points on the surface. The 
two-dimensional layer of atoms is then polycrystalline, however with 
almost certainly all of the crystals similarly oriented and representing only 
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lateral displacements with respect to each other without rotation. There 
can be only four, or possibly eight, different kinds of crystallites. Since 
all of these two-dimensional crystallites are simply related to the structure 
of the underlying nickel, and therefore to each other, it is not at present 
clear why they do not conspire to give a sharp pattern. 


Table 3. Patterns Representing the Q Structure 


(1) (2) (3) (4) (5) (6) (7) (8) 


Range Calculated beams 
Voltage 5 ae within, thie thee Observed beams 
and 
aes observ- ae 
length cee TEL. Spacing Azimuth+| Spacing Azimutht 
Fig. 9(A) 48 v 6A 
< 1-754 (20) 5-374 Ags 5-32 40 | A landsts 
(21) 4-06 4 19° 
(30) 3:59 A A;.B 
3-DA 
Fig. 9(B) 75 Vv 5-SA 
14154 (20) 5-374 A,B Just detectable 
(21) 4-064 Loe Just detectable 
(30) 3-59 A A,B 3-584 A, B weak 
(22) 3-104 C 3-104 C strong 
(31) 2-99 A 14° not seen 
PALIN 
Fig. 9(D) 120 v (10) 10-754 A,B ai 
1-124 (11) 6-22 4 C t 
(20) 5-374 ja gs) 5:59 A A,B 
(21) 4-064 19° 
(30) 3-59 4 A,B 3-704 A,B 
(22) 3-104 C 3-154 C 
(31) 2:99 & 14° 
(40) 2-69 4 A,B 2714 A,B 
(32) 2-474 Ze 2-474 23° 
(41) 2-354 PES 2-354 Ai ik 
(50) 2-154 A,B 2-154 A,B 
(33) 2-074 Cc 2-084 Cc 
(42) 2-03 4 EO® 2-054 20° 
(51) 1-934 9° 1-944 s° 
(60) 1:79 4 A,B I-824 A,B 
(43) 1774 25° 1-814 25° 
(52) 1-724 16° L774 16° 


} Azimuths other than that A, B and C are designated by the angle which they make 
with the A and B azimuths. 
t It should be pointed out that there are other diffraction patterns from the Q structure 


taken at 10 v and at 15 v (not reproduced here), which show the (10) and (11) reflections 
unmistakably and very strongly. 


With the crystal at room temperature, this type of structure seems to 
last indefinitely. If the crystal is warmed to 100°C above room tem- 
perature the diffraction spots become considerably sharper, presumably 
corresponding to crystal growth. “The spots can, however, never be made 
as sharp as one would expect them to be for the entire surface arranged in 


a single two-dimensional array. The patterns of fig. 7 are of intermediate _ 
sharpness. 
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At a temperature 225°c the patterns due to adsorbed gas atoms are very 
weak, and at slightly higher temperatures entirely missing, although the 
Laue spots from the nickel crystal are not greatly changed. This change 
with temperature seems to be due to melting the two-dimensional crystal 
of atoms so that they move about freely over the surface. A pattern 
reappears at once on allowing the crystal to cool. The pattern which 
appears after cooling is, however, never the 4 pattern. If the heating has 
been quite moderate, 225°, the new arrangement is that of the Q structure. 
If the heating has been more severe, at temperatures which may be as high 
as 400°C, the pattern may be the Q or the 3. There seems to be no heat 
treatment which can produce the reverse change, from the Q or 3 pattern 
back to the 4. Thus the arrangement of atoms represented by the 4 
pattern seems to be the least stable of the three arrangements. 

Whether the spots of a Q type diffraction pattern. are fuzzy, as in fig. 
9(C), or sharp, as in fig. 9(D), is determined by the temperature to which 
the crystal has been raised after deposition of the gas. Maximum tem- 
peratures of the order of 225°c lead to patterns with diffuse spots (like 
fig. 9(C)), but sharp spots are found after heating at 300°c or 400°c. 

Understanding the factors which determine the structure of adsorbed 
atoms is made difficult by inadequate independent control of temperature 
and gas pressure. The nickel crystal is heated by electron bombardment 
from a tungsten ribbon. Whenever this ribbon is heated the gas pressure 
rises (often by a factor greater than 10) because of the adsorbed nitrogen 
which it gives up, and whenever the ribbon is cooled there is a decrease of 
pressure (often by a factor of 10). A thin nickel crystal, which can be 
heated by conduction, is now being prepared. We expect that the factors 
which determine the different structures of adsorbed atoms can be found. 
out from experiments with this crystal. 


3.2. CO on Nickel 


Carbon monoxide molecules arrange themselves on the surface of the 
nickel crystal in the 4-structure. Except for this one point the behaviours 
of carbon monoxide and nitrogen are, however, quite different. The 
probability that a CO molecule which strikes the surface will stick to it is 
very high. Yet the CO is held very weakly, and molecules of CO are com- 
pletely removed by heating the crystal to only 175°c. Other differences 
in the behaviour of the two gases are the absence of polymorphic forms of 
two-dimensional CO crystals, and the greater tendency of CO molecules 
to arrange themslves in an almost perfect two-dimensional array, rather 
than imperfectly. A typical pattern from a monolayer of CO molecules 
is reproduced as fig. 10 (Pl. 30). . 

At an appropriate electron beam voltage and a suitably low CO pressure, 
the gradual development of a diffraction pattern can be observed quite 
clearly. An annular ring of nebulous character shows up first on the 
fluorescent screen. This ring seems at first circular, then assumes an 
angular appearance which soon can be identified as hexagonal. Finally, 
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the sides of the ring clear up and only six diffraction spots are left. Initially 
the spots are diffuse, but they rapidly become quite sharp. In the case of 
nitrogen, the development of a diffraction, pattern is less easily followed. 

The sticking probability measurements for CO, which are very crude, 
were made by observing the time required for the 4 pattern to appear at 
maximum intensity and sharpness upon the fluorescent screen immediately 
after the crystal had cooled from high temperature. A simultaneous 
record was made of the CO pressure, which was necessarily always changing 
rapidly. In a number of tests, values of the integrated product of the 
pressure and the time were found in the range 0-4 to 6 x 10-*mm Hg sec. 
The theoretical value for CO molecules striking the surface at 300°K 
and forming a complete monolayer with the 4 arrangement is 1-3 x 10~® 
mmHgsec for sticking probability unity. Although the observations 
were extremely crude, they do prove that the sticking probability is very 
high and of the order of unity. They indicate also that CO molecules are 
fairly free to move about over the surface at room temperature. 

The diffraction patterns from monolayers of CO molecules were always 
quite sharp and no change of sharpness was produced by warming the 
crystal after deposition or -by deposition upon a warm crystal. There 
was, however, no deposition of CO upon the crystal at temperatures in 
excess of 150°c, and this temperature was just high enough to clean the 
surface completely. 

In a series of interesting experiments, CO was deposited first upon a clean 
crystal surface, then the surface layer was removed by heating to 175°c 
after which a new layer of CO molecules was deposited. It was found 
that after two of these cycles the surface could be only partially cleaned at 
this low temperature, and the situation became worse after another cycle. 
The obvious conclusion from these tests is that we were observing the effects 
of impurities in the CO (from residual gas in the tube), which could not 
be removed at 175°c. Behaviour was normal again after the crystal had 
been heated above 800°c. 

In some experiments the Q structure was observed, but this was only 
after the crystal had been exposed over night to CO at an exceptionally 
high pressure (2 x 10-*mm Hg as an average value). It seems altogether 
probable that this structure was due not at all to CO but to impurities. 


3.3. Oxygen on Nickel 


The behaviour of oxygen is very different from that of nitrogen or CO. 
Exposure of the clean crystal surface to oxygen results in gradual weakening 
and final disappearance of Laue beams from the nickel, without the 
appearance of any new diffraction beams at low voltages. When a crystal 
which shows no diffraction patterns of any sort, after for example having 
been exposed to oxygen at the comparatively high pressure of 5 x 10-8 
mm Hg for 16 hours, is heated to progressively higher temperatures, the 
following series of changes can be observed. 
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After heating at 200°o, ditfuse low voltage beams can be seen (fig. 
11(A) Pl. 30) but the amount of oxygen on the surface is still so great that 
no Laue beams from the nickel are seen up to the limit of the observations 
at 400v. ‘Some sharpening of the oxygen diffraction pattern can be 
produced by further heating (fig. 11(B)), still without removing enough 
of the oxygen to permit the appearance of diffraction patterns from the 
nickel. More vigorous heating (probably+ above 500°C) results in re- 
appearance of the nickel diffraction beams at apparently full intensity, 
but with no change in the oxygen diffraction pattern which is now almost 
certainly the pattern from a single monolayer. Finally the diffraction 
pattern from the oxygen is destroyed by heating at a bright cherry red. 

The two dimensional spacings corresponding to the spots of fig. 11 are 
written down in table 4. The spots can be indexed as (22) and (40) reflec- 

tions of the Q structure (see table 3). We do not understand how a thick 
layer of annealed oxygen can give the same type of Q pattern which is 
produced by a single monolayer. 


Table 4. Patterns from Oxygen on Nickel 


Voltage Observed Beams Miller 
and —_ indices 
wave-length Spacing Azimuth (see table 3) 
Fig. 11(A) 82 Vv 2-66 A A,B (40) 
1-35 A 


3.4. Hydrogen on Nickel 


No evidence has been obtained of any diffraction pattern due to hydrogen. 
In one test hydrogen flowed through the experimental tube for 65 hours 
at an average pressure of 2x10-’mmHg. At the end of this time no 
low voltage diffraction features could be found; the Laue beams from the 
nickel lattice seemed to be even stronger than they were initially; and 
there was a considerable improvement in the apparent resolving power of 
the crystal, shown by each Laue beam being observable over a narrower 
range of wavelengths. It was shown, by observing a rise in pressure when 
the crystal was heated, that there was considerable adsorbed hydrogen. 

It seems altogether possible that hydrogen atoms or molecules may be 
as close together on the surface as are surface nickel atoms, and with the 
same lateral arrangement. If this were true, their presence would result 
only in modification of the intensities of diffraction beams attributed to 
the lattice of the crystal. The improved resolving power of the crystal 
remains to be explained. 


2 en 
+ The thermocouple attached to the crystal was out of order when this test 
was carried out. 


P.M. ; Q 
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§ 4. Discussion 


Many observations of relative intensities of diffraction beams from 
adsorbed gas layers require clarification. Most of this will have to await 
further work. 

Perbaps the simplest of these observations upon relative intensities is 
shown in the patterns of fig. 7. One notes that diffraction patterns repre- 
senting the 4 structure are strong at some voltages in the A azimuth and 
at other voltages in the B azimuth. “ It will be noticed that the nickel 
atoms in the surface are of two sorts—those that are adjacent to gas atoms 
and those that are not. One-fourth of the nickel atoms are of this latter 
class, and they together form a layer of the same structure, orientation 
and scale factor as the gas atoms. If the two classes of nickel atoms scatter 
radiation in different amounts—if, for example, the gas atoms shield the 
adjacent nickel atoms more effectively than the non-adjacent—then, of 
course, the layer of nickel atoms as a whole will give rise to a differential 
diffraction beam capable of interference with that due to the gas layer.’ 
This accounts for the variation of intensity with azimuth, and its change 
with voltage. In principle, one can work out from it the distance between 
the adsorbed nitrogen atoms and the surface nickel atoms. The electron 
wavelength is, however, too short to give any respectable precision to such 
a calculation. The data are consistent with a separation of 34 which is 
a reasonable value. 

It has already been pointed out that in different experiments, which 
have led to atoms arranged on the surface in the Q structure, relative 
intensities of diffraction spots have been widely different. This has been 
interpreted to mean that there are a number of different possible arrange- 
ments of atoms within the large two-dimensional hexagonal unit cell 
(dy = 12-4) of the Q structure. And even one single diffraction pattern 
must be interpreted to arise from a structure in which there are different 
arrangements within the unit cell on different parts of the effective surface. 
It is hoped and expected that the distributions of arrangements giving 
rise to particular patterns can be found out. 

In each diffraction, pattern from a monolayer of adsorbed atoms or 
molecules, of whatever sort, it has been observed in general that the pattern 
appears most strongly at some particular electron wavelength. A single 
diffraction spot of the pattern can always be traced over a range of electron 
wavelengths, satisfying always the appropriate plane grating formula, 
but its intensity varies greatly.. Such variation is not consistent with a 
pattern produced by a single atomically flat plane of atoms. The probable 
interpretation of this variation, which has not yet received theoretical 
quantitative confirmation, is that the variation is due to atomic steps on 
the surface within a distance appropriate to the wave-front of an electron. 
One single step, for example, would result in diffraction from two 
monolayers of atoms with a spacing of 2-03 A (the separation of (111) planes 


eT ee ne eM A 
} The quotation is from Davisson and Germer (1927). 
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of nickel atoms). This interpretation should be amenable to easy con- 
firmation. 

That this interpretation is correct is suggested by the experimentally 
observed fact that quite flat crystal surfaces are required for the observa- 
tion of any diffraction patterns from adsorbed monolayers. In earlier 
experiments, before the use of the present nickel crystal, tests were carried 
out upon a tungsten single crystal. Excellent Laue diffraction beams were 
obtained from this crystal. Single monolayers of adsorbed gas atoms were 
deposited upon its surface by well-known techniques, but in no case could 
any diffraction pattern at all be obtained from one of these monolayers. 
Tungsten crystals, of the type used in these experiments, are known to be 
notoriously rough. It is now believed that this was the cause of the failure 
of the tungsten experiments. It seems probable that the method of 
investigation outlined in this paper will be applicable only to selected crystal 
surfaces which can be made exceptionally flat on an atomic scale. 

The remarkable anomalous improvement in the resolving power of the 
nickel crystal, which is produced by hydrogen, is not understood. Is it 
possible that this improvement is the result of a change to normal spacing 
of an abnormal spacing between the first and second layers of nickel atoms? 
Accurate intensity measurements are required to settle this question. 


§ 5. New EXPERIMENTS 


A number of continuations of the present experiments immediately 
suggest themselves : 

1. The most pressing of these is the certain identification of the gas giving 
rise to the patterns assigned to nitrogen. This will require a nickel crystal 
which can be heated by conduction in order to avoid the fluctuations in 
pressure taking place on beating and cooling the tungsten ribbon used in 
bombarding the crystal. 

2. Such a crystal, which can be heated by conduction, should be valuable 
in other experiments also; for example, in understanding the factors 
which determine the transitions between the crystalline forms of the 4, 
Q and 3 structures, and by the flash filament technique in finding just when. 

we have a monolayer of adsorbed atoms in those cases in which it is now 
uncertain, and finally in determining unambiguously whether a gas is 
adsorbed as atoms or molecules. 

3. Repetition of some of these tests upon a (100) face of a nickel crystal. 
This would seem to require preparation of an unusual surface without 
(111) facets. Present electropolishing techniques offer promise that a 
satisfactory (100) surface can be obtained. 

4. Tests of the reaction of one gas with another upon a crystal ice. 
For example, will gas atoms (presumed to be nitrogen) adsorbed upon a 
nickel surface in the unstable 4 structure react more readily with hydrogen 
than will atoms in the other two more stable arrangements? Experiments 
of this sort require a considerable battery of valves with appropriate. 
auxiliary pumps and gas bulbs. These are now being installed. 


Q2 
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5. Studies of surface states of silicon crystals. The (111) surface of 
silicon can be made remarkably smooth on an atomic scale, and such 
surfaces should be very favourable also for studies of adsorption and 


catalytic reaction. 
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ABSTRACT 


A study has been made of the effect of 1 mev electron irradiation on the 
mechanical properties of lithium fluoride single crystals tested in compression. 
The yield point, c, has been. shown to vary with the irradiation dose, 4, 
according to the law c=o,,+ A[1—exp (—B¢)]}/2, where o,, 4 and B are 
constants. This relationship has been derived theoretically by considering 
the passage of dislocations through radiation induced obstacles dispersed in 
the slip plane, together with the capture of point defects by existing obstacles. 
The empirical relationship o=k¢1/?, previously used for radiation hardening 
in metals, is capable of representing only the medium dose results, but clearly 
does not fit at either low or high doses. : 

A new mechanism of cracking in heavily irradiated lithium fluoride, similar 
to that proposed by Cottrell for alpha iron, has been observed. With com- 
pression along a [001] axis, cracks are initiated at the intersections of 
orthogonal {110} type slip planes and extend on the (100) cleavage plane 
passing through the intersection and parallel to the compression axis. 


§ 1. InTRODUCTION 


Two types of irradiation hardening are theoretically possible in crystalline 
solids which normally deform by a slip process. Either (a) the nucleation 
of slip may be made more difficult by making dislocation sources harder 
to operate (source hardening), or (b) the stress for the propagation of slip 
may be increased (lattice hardening) by an increase in the stress required 
to move dislocations through the crystal. In well-annealed metals such 
as copper and nickel both types of hardening occur after neutron irradia- 
tion. In lithium fluoride the existence of lattice hardening after electron 
irradiation has been proved by observing the dislocation patterns (revealed 
by etching) produced by hardness indentations (Whapham 1958). In 
unirradiated material dislocations glide appreciable distances away from 
the indentation whereas after irradiation dislocations can be seen only in 
the immediate vicinity of the indentation, the irradiation induced lattice 
hardening preventing these dislocations from moving away. 

Work preliminary to the present investigation has shown that it is the 
stress necessary to move dislocations that governs the plastic deforma- 
tion of lithium fluoride. The mechanism of deformation of stressed 
crystals was investigated by etching to reveal the slip dislocations. It 
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was found that in all specimens, both unirradiated and electron-irradiated, 
individual slip planes extended from their point of origin slowly, with 
progressive increase of stress. If the operation of dislocation sources 
determined the stress at which plastic deformation, occurred then complete 
slip lines would be expected to form abruptly. This work also showed 
that dislocations in the compression specimens move at stresses well below 
the macroscopic yield stress at normal strain rates. In unirradiated 
material dislocation movement is first seen at a resolved shear stress of 
0-13 kg/mm?, whereas the macroscopic yield occurs at a shear stress of 
0-5 kg/mm? when the applied strain rate is 10-*persec. In the electron- 
irradiated material, dislocation movement also starts at stresses roughly 
0-4kg/mm? less than the value at the macroscopic yield. 

Gilman and Johnston, who have measured the yield stress of lithium 
fluoride during bend tests (Gilman and Johnston 1956, 1958), and the 
stress dependence of dislocation, velocity as a function of neutron irradia- 
tion (Johnston and Gilman 1958), concluded that the irradiation hardening 
is governed by the resistance to motion of dislocations through the crystal 
lattice. In addition the form of the stress-strain curve can be explained 
in terms of the variation of dislocation velocity with stress and the multi- 
plication of dislocations. It is therefore apparent that the yield point at 
a given strain rate is a measure of the lattice hardening. No direct 
measurement of an increase in stress required to operate dislocation 
sources in lithium fluoride has yet been made. 

The present experiments were undertaken to determine in detail the 
dependence of the lattice hardening on the irradiation dose, and to com- 
pare this with a theoretically predicted relationship. 


1.1. Dependence of Lattice Hardening on Irradiation Dose 


Lattice hardening must be caused by ‘obstacles’ of some kind distri- 
buted throughout the crystal lattice which oppose the motion of dis- 
locations across the glide plane. It is not necessary to specify at this 
stage the nature of these obstacles except to point out that they are 
unlikely to be individual point defects during room temperature irradiation 
in view of the low temperatures at which these become mobile in most 
crystalline materials. It is more probable that each obstacle is a cluster 
or aggregate of defects of some kind. 

If N is the number of dislocations per unit volume held up against V, 
obstacles per unit area of glide plane then the strain rate € ata tempera- 
ture 7’ and stress o is given by 


. Nbv U(o 
é= Syetom(— a). i Ser PT ae 


where U is the activation energy for cutting through the obstacles, b is 
the Burgers vector and v» the frequency with which a dislocation line 
attempts to pass the obstacle. Solving eqn. (1) for o gives the lattice 
hardening stress as a function of Nz. To do this is is first necessary to 
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determine the stress dependence of the function U(c). Various equations 
can be chosen to represent the energy profile of the obstacle which is 
expected to be as shown in fig. 1. Application of stress to the crystal 
modifies the shape of the potential energy—distance curve as in fig. 2 by 
introducing a variation in potential energy with distance of ol,'b where 
I, is the average distance between obstacles in the slip plane. Mott and 


Fig. 1 
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DISTANCE THROUGH OBSTACLE. 


Energy profile of a dislocation cutting through an obstacle at zero externally 
applied stress. 


Fig. 2 
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‘Energy profile of a dislocation cutting through an obstacle when an external 
stress a is applied. 
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Nabarro (1948) and also Seeger (1958) have shown that a good approxi- 
mation to the new value of U is given by 


o \382 
v=0,(1-2) 


where U, is the activation energy at zero stress and oy is the stress required 
to overcome the barrier without thermal activation. 


When U=0, 
, Uy 
2x,bl,o\32 
and hence U=U,(1- ee) 23>. ne 
0 


Seeger derives an almost identical expression for U except for the constant 
4 instead of 2. 


In the simplest case, where no allowance is made for a stress dependence 
of 1, 


j 1 
ly = N22 
2Qxyba \32 


Substituting for U in eqn. (1) 


U kT Nb, \ ) 28 
— 12 o 1 = ae bed ae 
cuit leering) al 


is obtained and if NV, is proportional to the irradiation dose, ¢, initially, 
oo pi?, 
Friedel (1956) has suggested that J)’ should decrease with increasing stress 


due to the greater probability of a dislocation bowing out sufficiently to 
meet another obstacle at high stresses. Friedel calculates that 


Inserting this into eqn. (2), we obtain 


2: 4/3 1/38,42/8 3/2 
he u.|1- = zoe =r | 


U Nz 1/3 


Substituting this modified expression for U into eqn. 1, the expression for o 


becomes 
Nz 3/2 Nbv,\) 22732 
‘lel lea! Dt Gey 


In this expression also 
oa N12 


In both cases therefore dislocation theory predicts a dependence of lattice 
hardening on Nz". 
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Initially the number of obstacles NV, is expected to be proportional to 
the irradiation dose 4 but, as the concentration of obstacles increases, 
there is an increasing probability that new defects will be captured by 
existing obstacles rather than nucleating new ones. Hence the number of 
new obstacles formed during a given increment of dose decreases as the 
dose increases. If there is a volume v associated with each obstacle within 
which capture will occur, then the rate at which new obstacles are formed 
N,, can be written 

N,=ad¢(1—vN,), 


where « is the number of obstacles produced per electron and ¢ is the dose 
rate. Integrating, the number of obstacles N, is given by 


= 1 
N= - [1 —exp (—avd)] 
and hence the lattice hardening is given by 
if 
o;C 5 [1 —exp (—av¢d) }!? 
or o;=A[1—exp (— B¢)}?”. BE os STAG) 


This formula predicts that initially the lattice hardening o; will be 
proportional to #1? and that at large doses o; will saturate at a constant 
value A. The relationship does not depend upon the properties of any 
particular material and is intended to apply generally to crystalline solids 
which deform by slip. 

The alternative relationship which has been fitted to the hardening of 
copper single crystals by neutron, irradiation (Blewitt 1956) is that 

a,x pt, 

This relationship has no sound basis in dislocation, theory and is purely an 
empirical one which happens to fit the experimental data over the range 
of doses used. The importance of deciding between, these two expressions 
is of considerable technological as well as fundamental interest since the 
¢'3 relationship predicts a continually increasing hardening with neutron 
exposure instead of the saturation expected in the new theory. In metals 
very high neutron doses are required to decide between the two theories 
since at low doses by a suitable choice of constants the two relations are 
very similar. In lithium fluoride damaged by electron irradiation however, 
only comparatively short irradiation times are required to introduce 
enough damage to decide which relationship is more correct. 


§ 2. EXPERIMENTAL 


Specimens, 2-4 x 24x 20mm in size and bounded by (100) cleavage 
surfaces, were cleaved from two single crystal blocks 1$in. x lyin. x fin. 
selected from a batch of material supplied by the Harshaw Chemical 
Company. The cleaved specimens contain fresh surface dislocations 
caused by the cleavage process and dust scratches, together with about 10° 
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aged dislocations per cm?. Roughly half the aged dislocations are con- 
tained in low-angle boundaries and half are randomly distributed through 
the crystal. The two crystal blocks used in this work were of similar 
purity and had been subjected to similar thermal treatment and the yield 
stresses, which are very sensitive to these factors, were practically 
identical (0-498 + 0-010 and 0-493 + 0-010 kg mm? respectively). 

The specimens were irradiated, three at a time, in a beam of 1 Mev 
electrons produced by a Van de Graaff generator. As the penetration of 
the electron, beam in lithium fluoride is only 1-2 mm and specimens thicker 
than this were required in order to perform satisfactory tests, the 2-4mm 
thick specimens were irradiated firstly on one face and then turned round 
to irradiate the opposite face. Using this technique the radiation dose 
was fairly uniform throughout the crystal, as shown by optical absorption, 
and consistent results were obtained with them. Crystals thicker than 
those used show a light band in the centre, where the irradiated regions 
do not overlap and thinner ones show enhanced absorption in the centre. 

The maximum dose used was 101° electrons/cm? since, above this dose 
cracking, produced by radiation -expansion of the lattice, was encountered. 
It is unlikely that strains present in specimens which had not cracked 
could affect the yield point since most of the stress-strain curves showed 
an easy glide region of 14% even at the highest doses used. This deforma- 
tion should eliminate any non-homogeneous strain present. 

Care was taken to ensure that the temperature rise during irradiation 
was not excessive. Variation of the beam current below the maximum of 
0-54 .A/em? did not have any observable effect in subsequent testing. 
Theoretically it is estimated that the maximum temperature rise was 
about 10°c. 

After irradiation, steel end pieces and two brass cross pieces, 10mm 
apart, were fitted to the specimens, using a synthetic resin adhesive. 
Mounted in this way they were placed in a special compression jig designed 
to produce axial loading and tested at a strain rate of 10-*sec. The strain 
of the 10mm central section of the specimens was measured by a capaci- 
tative strain measuring head attached to the two brass cross pieces. With 
this arrangement strains from 0-01°% to 10% could be measured. 

Representative stress-strain curves are shown in fig. 3. The 
unirradiated crystal yields at 0-496 + 0-010 kg/mm2. Out of ten tests, three 
of the stress-strain curves show a flat easy glide of 0-3, six showed a 
rate of hardening in this region which is between one-quarter and three- 
quarters of the following linear hardening rate, while one did not show any 
sign of an easy glide region. The rate of work hardening in the linear 
hardening region varied from test to test between the limits 0-07 and 
0-18kg/mm? % strain. The curves with the flat easy glide region are very 
similar to that shown by Gilman and Johnston, (1956). 

After 10% strain the specimens were considerably bulged. Though, in 
the early stages of deformation, slip on all systems could be seen (by 
etching to reveal dislocations), in heavily strained specimens gross slip 


Hardening of Lithium Fluoride by Electron Irradiation 243 


occurred only on one or two systems since, while one pair of crystal faces 


ee considerably bulged, the pair at right angles were macroscopically 
at. 
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Stress—strain curves of irradiated lithium fluoride specimens tested in 
compression. 


It is seen from the graphs that 1 Mev electron irradiation considerably 
increases the yield point. A dose of 10” electrons/cm? just produces an 
‘increase clearly beyond experimental error and at 10'*electrons/cm? the 
yield stress is raised to 15 times the unirradiated value. The easy glide 
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region increases to 14% strain, even at a dose of 10™ electrons/cm? and 
remains at approximately this value for doses up to 3 x 1014 electrons/ cm? 
though, as in the unirradiated case, the occasional test will show no sign of 
easy glide (e.g. the test illustrated in fig. 3 for a specimen irradiation of 
1013 electrons/em3). At doses above 3 x 101* easy glide regions of several 
per cent were obtained (8% in the case of one specimen which had been 
given a dose of 10" electrons/em?). One specimen which had been given 
a dose of 6:25 x10 electrons/em?, showed an easy. glide region of 
about 4% and in the subsequent stage of up to 10% strain the hardening 
was only 0-05kg/mm2 % strain, whereas another specimen did not show 
any easy glide but a hardening rate of 0:15kg mm? % strain. Onexamina- 
tion of these two specimens after testing, the one which gave 4% easy 
glide showed deformation predominantly on only one set of {110} slip planes 
whereas the other specimen, showed considerable intersecting slip. This 
dependence of the hardening rate on the degree of intersecting slip is well 
known in unirradiated materials. On deformation of irradiated crystals 
fewer slip lines are formed for a given strain but each line contains a higher 
linear density of dislocations than is found in the unirradiated crystals. 
After irradiation, therefore, a greater strain is expected before slip on 
intersecting planes produces hardening. This is in agreement with the 
longer easy glide region in irradiated crystals. 

No cracks were observed in unirradiated crystals and crystals given an. 
irradiation dose of up to 3 x 10! electrons/cm? when they were tested up 
to 11% strain. Cracks were first observed at 10$°% strain in specimens 
given, a dose of 1018 electrons/cm?, 6°% strain with a dose of 3 x 104, 3% 
strain at 1-5 101° and a specimen irradiated to 3x 10% electrons/em? 
which had been given only 0-3°% strain showed many cracks. Even in 
these highly irradiated specimens which developed cracks in the early 
stages of deformation, strains of 1 or 2°4 could be obtained before the 
compression specimens shattered. The mechanism of cracking in the low 
dose crystals (< 101° electrons/cm?) is complicated, the cracks not being 
confined to the {001} cleavage plane. After high doses, however, a new 
mechanism appears to operate. The cracks occur entirely in the {001} 
type planes and are initiated at the intersection of orthogonal slip planes as 
shown in figs. 4.and 5 (Pl.31). Even at a magnification of 500 there is no 
evidence of cracking on other planes at the point of initiation. The cracks 
occur in the crystal in the form of slits, which sometimes extend to the 
crystal surfaces. This type of cracking is different from that observed by 
Stokes e¢ al. (1958, 1959) in unirradiated MgO and is similar to that 
proposed by Cottrell (1958) for alpha iron. In this mechanism slip 
dislocations on orthogonal {110} type planes coalesce to form a new 
dislocation on the intersection of the planes. In fig. 6 the dislocation 
reaction is 


< [101] + < (101]->a[001]. 


The new dislocation lies in the (100) cleavage plane and is equivalent to an 
extra half plane inserted between (100) planes. Other slip dislocations 
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following the original ones may also react together to form more a[O01 | 
dislocations so increasing the number of extra half planes on one side of the 
intersection until a crack forms as in fig. 6. 


Fig. 6 


loo) 


jae ~OMPRESSION -DIRECTION—— 


Coalescence of slip dislocations on two slip systems to form new dislocations 
with extra half planes on the (100) cleavage plane to form a (100) crack. 


An interesting feature of crystals which have been irradiated to doses of 
10" electrons/cm? and above is that the slip traces show an enhanced optical 
absorption. This is seen clearly in fig. 4. The enhanced darkening of 
deformed crystals on irradiation has been reported, but this is the first 
observation of enhanced darkening by deformation of irradiated crystals. 
If, for example, the irradiation dose given to the crystal has produced a light 
brown, colouration in it, slip produced by subsequent deformation is 
observed as dark brown planes extending through the crystal. The 
optical density in other regions, even near the slip planes, is not changed. 
It is possible that this enhanced optical absorption is due to colloidal 
lithium. 

The rate of work-hardening in the linear hardening stage of the stress— 
strain curves for irradiated erystals is essentially the same as in the 
unirradiated case. There is a considerable scatter of individual hardening 
rates however, which would obscure a small effect. 

The mean of the critical resolved shear stresses for the easy glide region, 
is plotted against the irradiation dose ¢ in fig. 7, against ¢'” in fig. 8 and 
é'3 in fig. 9. Only in fig. 8 is the initial part of the curve linear indicating 
that at low doses (< 10! electrons/em?) the correct relationship for the 
dose dependence of the critical shear stress is 


o oc pi, 


At doses greater than, 10' electrons/cm?, however, the critical shear 
stress increases less rapidly than the 4"? relationship and finally reaches 
a saturation value of 7:-63kg/mm?. After 5 x 101° electrons/cm? the critical 
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Plot of the critical resolved shear stress versus irradiation dose. The relation 
o=0-5+7-13[1—exp (—6-18 x 10-164) 2 is shown as a solid line; the 
broken line is g=5-18 x 10-541/3, 
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Plot of the critical resolved shear stress versus the square root of the irradiation 
dose. 
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shear stress is within experimental error of the saturated value. The 
results have been fitted to an equation of the type 


o=0,+0,=0,+A[1—exp (—B¢)|'? 


where o,,=critical shear stress of unirradiated material and oc; has been 
derived theoretically earlier in the paper. 
This gives 


o=0:5+ 7:13[1—exp (—6-18 x 10716) }¥? ko/mm?. 


This relationship is plotted as a continuous line in each graph and rep- 
resents the experimental points very accurately. 
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Plot of the critical resolved shear stress versus the cube root of the irradiation 
dose. 


It is seen from fig. 9 that the previously used relation to describe the dose 
dependence of irradiation hardening in copper, namely 


ox fis, 


does represent the experimental results between ¢ = 10" electrons/cm? and 
3x 10% electrons/em?. An, incubation period is observed at low doses 
however, and above 4 x 10! electrons/cm? it becomes quite clear that the 
cube root relationship is not the correct one. 
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§ 3. Discussion 


The most important result of this work is the establishment of the 
equation o=0, + Af 1—exp (—B¢)}? 


as the correct relationship between irradiation hardening and electron 
dose in, lithium fluoride crystals. This equation has been derived from 
dislocation theory by considering the passage of dislocations through 
obstacles dispersed in, the slip plane coupled with a saturation in the 
number of obstacles due to the probability of the capture of point defects 
by existing obstacles. At present it is not possible to determine unequi- 
vocably the volume v over which capture will occur because in the equation 
only the product av can be measured. An estimate of the value of a 
depends upon the number of defects which form an obstacle and the 
number of defects produced per unit irradiation dose. 


Number of Number Volume of Radius of 
. oint of sphere of sphere of 
Ae e ee clusters influence influence 
ee Le ina per around a around a 
Tea: cluster electron cluster, cluster 
(1%) (a) (vin 10-18em3) (7) 
l 200 3 9OA 
2 100 6 L1OA 
200 5 40 15 1604 
10 20 30 1904 
20 10 60 2404 


Table showing how the volume v and radius r of the sphere of influence 
around a cluster varies with the value assumed for the number of point defects 
in a cluster n,. 


Values of » and the corresponding radius 7 based upon, obstacles of 
various numbers of defects are given in the table. In this table it has been 
assumed that the cross section for the production of point defects in lithium 
fluoride is the same as the value (3000 barns) measured for the production 
of F centres in KCl by electron irradiation (Varley 1954a, b, 1956). The 
number of point defects estimated in this way will be a lower limit as it 
is possible that there may be some negative ion vacancies without an 
associated electron, in addition to the F centres. Here n, 18 the minimum 
number of point defects which can form an obstacle which will oppose 
dislocation movement. It is an essential part of the theory that new 
defects will attach themselves to these obstacles making them grow in 
size as the irradiation dose increases. The volume v will be constant 
throughout the irradiation and will correspond to the minimum size 
obstacle. It is seen from the table that the most probable value of r lies 
between 100 and 250 4. 
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It is not certain whether individual point defects could produce hardening 
(n,=1). Gilman and Johnston (1958) notes the correlation between the 
onset of radiation hardening and the appearance of the 4500 4 absorption 
band in neutron irradiated lithium fluoride. This 45004 band has been 
ascribed to M centres (F-centre +metal-halogen vacancy pair) or some 
sort of electron excess centres which are larger than F centres. The 
etching effect in both neutron and electron-irradiated lithium fluoride 
observed by Gilman and Johnston (1958) is strong evidence that the 
obstacles grow in size with irradiation. A threshold dose is necessary tor 
this etching effect to occur. This would be expected if obstacles could 
only produce an etching effect after they had grown to a critical size. 

The investigations of the small angle scattering of x-rays by neutron 
irradiated lithium fluoride by Lambert and Guinier (1957 a) and Smallman 
and Willis (1957) show that small cavities having an average radius of 
gyration of 10-20A are present after exposures of greater than 5 x 1016 
neutrons/em?. Though no small angle scattering work has been done 
on electron irradiated lithium fluoride the larger cavities produced by 
subsequent annealing were observed and are shown, in fig. 10 (Pl. 32). 
Further evidence of the nature of the defects in irradiated lithium fluoride 
is provided by the observation, by Lambert and Guinier (1957 b) of diffuse 
x-ray streaks along {100} directions in neutron irradiated material. These 
authors interpret these streaks to be caused by thin platelets of lithium 
metal lying on {100} planes. 

The lattice hardening introduced by irradiation is on the new theory 
additive to the initial critical shear stress of the material. Careful analysis 
of the results show that the experimental formula is of the form 


o=o,+f(¢) andnot o=f(¢+¢,), 


where ¢,, is a ‘ virtual’ dose required to produce the critical shear stress of 
the unirradiated crystal. This indicates that lattice hardening introduced 
by irradiation is fundamentally different from the hardening responsible 
for the initial critical shear stress of the material. 

The new relationship derived for the dose dependence of the lattice 
hardening component of irradiation hardening, 


o;= A[1—exp (— Bd) }!” 


has been verified experimentally only for lithium fluoride at present. It 
is expected that the formula will apply generally to crystalline solids which 
deform by slip, however. Preliminary results on, single and polycrystals 
of copper indicate that the equation fits the results up to the highest doses 
so far tested (102° neutrons/cm?). The divergence between eqn. (3) 
and a,x ¢¥3 becomes large, however, only at high doses and the present 
data does not yet enable a definite decision to be made between the two 
relationships. 
The decision between these theories is of considerable technological 
importance since the 4"? relationship predicts a continually increasing 


PM. E 


250 On the Hardening of Lithium Fluoride by Electron Irradiation 


hardening with dose instead of saturation expected in the new theory. 
The two theories make considerably different predictions, therefore, when 


extrapolated to very high doses. 
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ABSTRACT 


The scattering of electrons and x-rays by stacking faults is due to a phase 
change in the diffraction by erystal planes which do not lie parallel to the 
fault. A simple model of the process in copper estimates a resistivity 
p=1x10-'? Bohm-cem for a total stacking fault density of Bem on all 
planes. A wave matching calculation in which Tamm surface waves are used 
confirms the order of magnitude of the effect and is also in agreement with the 
prediction of the simple model that the distance over which scattering from 
two faults could interfere is less than the lattice spacing. These results are 
consistent with the available experimental data and it is concluded that the 
effect explains the resistivity caused by extended dislocations in cold worked 
metals and by tetrahedra of stacking faults in quenched metals. 


§ 1. [INTRODUCTION 
1.1. Resistivity of Dislocations : 

Experimental Data and Theoretical Predictions 
DisLocaTions and point defects created during plastic deformation of 
a metal cause an increase in the residual electrical resistivity. The part 
of this resistivity increase which remains after annealing at temperatures 
at which point defects can diffuse is ascribed to the scattering of conduction 
electrons by dislocations. 

From the results of Clarebrough et al. (1955) for the changes in density, 
stored energy and electrical resistivity of cold-worked copper during 
recrystallization, Seeger and Stehle (1956) quote for the observed resis- 
tivity increase a figure p=2-5x10-*ohm-cm per 1% density change. 
If we assume that each atomic length of dislocation in the crystal increases 
the volume by 1-5 atomic volumes (Stehle and Seeger 1956) the dislocation 
density to be associated with this change in resistivity is approximately 
1013em-2. Thus the experimental figure for the resistivity increase in 
copper due to a dislocation density of N em is p= 2-5 x 10~! N ohm-cm, 
More recent results of Clarebrough et al. (1957) for copper deformed in 
compression, confirm this result and also indicate that the dislocation 
densities deduced by the above method from the density change are in 
reasonable agreement with those deduced from the stored energy release 


on recrystallization. 
ne EE eee 
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The scattering of conduction electrons in copper by dislocation strain fields 
has been calculated by Hunter and Nabarro (1953) and for a dislocation 
density of N cm they find a resistivity increase p= 3-8 x 10~*! V ohm-cm 
which is almost two orders of magnitude smaller than the experimental 
figure given above. A slightly higher estmate of p= 6-0 x 10°! V ohm-em 
was obtained by Seeger and Stehle (1956) by considering the additional 
scattering effects due to anharmonic displacements and the resulting 
redistribution of electric charge. There might also be considerable 
scattering of electrons by the dislocation core where the atomic positions 
and electron charge distribution are not well known. Harrison (1958) has 
calculated the scattering from a hollow core and obtains a resistivity 
p=5x 10-79 N ohm-cm. 

It appears therefore that even with such an extreme model of the core 
it is still not possible to account for almost four-fifths of the experimental 
value for the resistivity. Later we shall quote an example where scattering 
from the dislocation core is even less able to explain the observed results. 


1.2. Resistivity of Stacking Faults: Experimental Evidence 


A possible explanation of these anomalous dislocation resistivity results 
was first advanced by Broom (1952). Dislocations in the face-centred 
cubic lattice may dissociate, in a manner proposed by Heidenreich and 
Shockley (1948), into two partial dislocations separated by a ribbon of 
stacking fault whose width is determined by the balance between the 
repulsion between the two partials and the surface tension of the stacking 
fault. Broom suggested that stacking faults may be effective in scattering 
electrons and that in the case of wide ribbons most of the observed resis- 
tivity might be due to the faults. A study of the increased resistivity of 
nickel-cobalt wires of various alloy compositions after a standard cold 
drawing treatment (Broom and Barrett 1953) showed a maximum resistivity 
for the alloy near the transformation composition where the stacking 
fault energy might be expected to be low and the fault ribbons wide. 
If we can asume that the dislocation and point defect densities in all the 
wires were similar this is strong evidence in favour of a large resistivity 
due to scattering from stacking faults. Broom and Barrett were also able 
to confirm the presence of wide faults from x-ray diffraction measurements. 
The displacement and broadening of the diffraction peaks as a result of 
the presence of stacking faults has been calculated for a random distribution 
of parallel faults by Paterson, (1952) in terms of a faulting probability «. 
Unfortunately the interpretation of « is obscure when the faults have a 
limited width. Nevertheless it has proved possible (Christian and 
Spreadborough 1956) to establish a strong correlation in several alloys 
between the manifestation of a large resistivity increase after cold work and 
a large probability of faulting « in alloy filings. The experimental evidence 
in favour of a large resistivity due to stacking faults produced in plastic 
_ deformation is thus fairly convincing if not altogether conclusive. 
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To obtain an estimate of the stacking fault resistivity required to 
explain the experimental results we shall assume that scattering due to 
the dislocation core and strain field is correctly given by the theoretical 
predictions quoted above and is thus responsible for not more than one-fifth 
of the total effect. If we assume the average stacking fault ribbon width 
in copper (Seeger and Schoeck 1953) to be 254 it follows that with a 
dislocation density of V cm-* there is associated a stacking fault density 
B=2-5x 10-7 Nem considered to be equally distributed over all four 
slip planes. Substituting for NV in the experimental value for p given 
above we find for the resistivity due to such a fault density, 8 (when the 
resistivity is entirely due to faults) 


p=1x10-"Bohm-em. At wemo cs ran gerne 


This figure is however subject to serious qualifications in view of recent 
suggestions (Seeger ef al. 1959) that the stacking fault energy of copper 
is much higher than hitherto supposed. The equilibrium width of the 
stacking fault ribbon would then be much less than 254. Seeger (private 
communication) has also made new calculations of the resistivity due to 
the dislocation strain field and finds that the estimate should be considerably 
increased. The effect of these two changes is apparently to leave the 
predicted value for the resistivity due to stacking faults unchanged. 
In view of these uncertainties it is fortunate that estimates of the stacking 
fault resistivity can also be obtained from observations on, quenched gold. 

A critical study was made by Bauerle and Koehler (1957) of the annealing 
of the resistivity increase observed in gold wires after rapid quenching 
from high temperatures. It seems probable that the part of the resistivity 
which anneals out at temperatures above 500°C is due to scattering from 
tetrahedra of stacking faults of the type observed in foils of quenched 
gold by Silcox and Hirsch (1959) using the transmission electron microscope 
technique. Sileox (private communication) finds, on the basis of Bauerle 
and Koehler’s data and the microscope observations, that a residual 
resistivity of p=6 x 10-°ohmem, found after quenching from 900°C and 
then annealing at 40°c, is to be associated with a total stacking fault 
density B=1-:2x104cm- arranged in the form of tetrahedra. The 
resistivity of stacking faults is then 

GU PoqlOs AR OHM-cM es faab fee \anhs ee a (4) 
where f is again the total fault density on all four slip planes. 

It shouldbe noted that almost all the resistivity due to the tetrahedra 
of stacking faults must come from the stacking faults themselves. The 
dislocation line length (along the edges of the tetrahedron of average 
length 300A) associated with a given area of fault is smaller by a factor 
of almost ten than the dislocation line length associated with the same 
area of fault in the deformed metal. In addition the dislocations on the 
edges of the tetrahedron have the small Burger’s vector ¢ [110] and are 
arranged so that their long range stress fields cancel. The conclusion is 
that the dislocations associated with the tetrahedra cannot contribute 


254 A. Howie on the 


more than a negligible amount to the scattering. Therefore although 
estimate (1) for the stacking fault resistivity is of more direct interest since 
it refers to copper (the metal to which the calculations given later apply) 
the quenching data for gold are probably the most reliable evidence in 
favour of a stacking fault resistivity of this order of magnitude. 


1.3. Theory of Resistivity of Stacking Faults 


A critical survey of previous attempts to investigate the scattering of 
conduction, electrons by stacking faults has been given by Seeger (1956). 
He points out that in the free electron approximation (in which the 
potential is assumed to be constant inside the metal) electrons do not feel 
the presence of a stacking fault at all. For this reason a treatment by 
Klemens (1953), in which the fault was replaced by a scattering potential 
and a large reflection coefficient R=4 estimated for electrons incident 
normal to the fault, cannot be valid since free electrons are scattered by 
such a potential. Another estimate of R=0-4 by Klemens (1956) was 
obtained by means of a deformation potential method but this approach 
holds only for slowly varying strains and is far outside the range of its 
applicability when used in the case of a stacking fault where the shear 
is infinite at the fault and zero everywhere else. An alternative approach 
to the problem is to use Bloch wave functions in the perfect crystal on 
either side of the fault and match these functions at the fault. Such 
methods have been used by Blatt et al. (1956) and Tweedale (see contri- 
butions of Ziman, to the discussion at the end of the papers by Klemens 
(1956) and Seeger (1956)). Tweedale was able to perform a more complete 
calculation and estimated a very small reflection coefficient R of the order 
of 10-3 to 10-4. However, as Seeger has pointed out, in neither of these 
investigations were surface waves at the fault considered and these waves 
are an essential feature of the problem. 

Seeger considered the scattering from the fault in terms of a diffraction 
mechanism and obtained a reflection coefficient R between 4 and 4. 
Unfortunately the calculation entailed a rather implausible use of pertur- 
bation theory (see § 3.1) and certain important principles of the diffraction 
process are obscured. 

Ziman, (1958) has also emphasized the part played by diffraction in 
the problem and has predicted a reflection coefficient proportional to the 
intensity of the diffracted crystal wave. He has pointed out that in view 
of recent results in, connection with the Fermi surface in copper (see § 2) 
the intensity of the diffracted beam may be large. However, he has not 
considered in detail the effectiveness of the various diffracted beams in 
contributing to the resistivity. 

In this paper the scattering from stacking faults as a diffraction 
phenomenon is examined more closely. In § 2 the necessary preliminary 
description of the part played by diffraction in the propagation of conduction 
electrons is given in terms of a specific model of the Fermi surface in, copper. 
It then proves possible (see §3) to build up a simple physical picture of 
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the effect of the fault and, by considering the equilibrium form of the 
crystal waves at large distances on either side of the fault, to derive an 
expression for the scattering and hence to deduce the resistivity. In § 4 
a wave matching treatment is employed using both Bloch waves and 
surface waves and an estimate for the transmission coefficient of the 
fault is obtained. The results of these calculations are discussed. 


§ 2. THe Fermi Surrace 1n Copper ; Two Beam APPROXIMATION 


It will be necessary in the calculations which follow to have a description. 
of the conduction electrons in copper. The Fermi surface in copper has 
recently been determined by Pippard (1957) on the basis of the anomalous 
skin effect. Figure 2 shows a section and fig. 3 a stereographic contour 
map of the surface which is approximately spherical, as would be expected 
in a monovalent metal, but is attracted to the nearest faces (the (111) 
faces) of the Brillouin zone rather more strongly than had been generally 
supposed. This means that for conduction, electrons travelling near the 
[111] direction diffraction from (111) planes in the crystal is important. 

It is proposed in the following treatment to apply the approximation 
of only two strong crystal beams—a direct beam and a diffracted beam 
from one set of crystal planes. The wave function is then written 


%(r) =exp (vk. r) + C(k) exp [7(k +g). r] ee eS) 


where k is the wave vector |k|=27/A, A being the wavelength, g/27 is a 
reciprocal lattice vector (in our case g=2z/a [111] and C(k) is the ratio 
of the amplitudes of the diffracted and direct beams and thus a measure 
of the importance of diffraction (or the effect of the Brillouin zone). If the 
crystal potential is given, by 


: h? : 
Vir)=>, V, exp (ig.r)= Boe. U, exp (2g. r) ee (4) 
the energy Z, inside the Brillouin zone, is obtained from the expression 
2 
= aS —U, 
= 
= ${k? + (k+g)}— f{Le{(k+gP—-k3 P+ [UF  - - ©) 
and the constant C(k) in eqn. (3) can be written in the alternative forms 
k2—y 
O(kj=—- ——, (6) 
l = SE 
U. 
C(k) = — a (7) 


(k+g)?—k?} + V/{[H(k+g)?—k?} P+ [U,P} 
These two expressions for O(k) will only be equivalent when the two-beam 
approximation is valid and all other diffracted beams have negligibly 
small amplitudes. It should be emphasized that this condition does not 
necessarily correspond to almost free electrons. In fact the two-beam 
approximation is best in the vicinity of the Brillouin zone where the 
amplitudes of the two beams are nearly equal. 
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We shall assume that the surface near the zone has the form deduced 
by Pippard using a one dimensional Kronig~Penney model and shall use 
eqns. (6) and (7) to obtain values of O(k) for wave vectors k on the Fermi 
surface where (f2/2m)y=5-56ev. At the Brillouin zone |C(k)|=1. This 
follows automatically from eqn. (7) but not from eqn. (6) unless |U,| is 
chosen correctly. From the value of k where the surface meets the zone 
we find (A?/2m)|U,|=3-25ev so that the energy gap at the zone is 6-5ev 
and compares well with Pippard’s prediction of 7-5ev on, the basis of the 
Kronig—-Penney model. We can now write eqn. (7) in a more useful form 
by expressing d, the distance from the end of the vector k to the Brillouin 
zone measured normal to the zone, as a fraction € of the distance [L = 4|g| 
(fig. 2). Then (k+g)?—k?=2|g|d= eg? and using the figure deduced 
for (h?/2m)|U,| and also the value for (h?/2m)g?=34-4ev eqn. (7) can be 
written 

1 
ON be (5-3) + 4/[(5-3e)? +1] (8) 

The surface can then be divided into rings (fig. 2) which are parallel 
to the zone and correspond to increments of «. A value of |C(k)| can then 
be obtained for the middle of each ring using eqn. (8) and these are the 
first of the three values of |C(k)| for each ring given in table 1. The second 
and third figures for each ring are obtained from eqn. (6), the second from 
points k on the surface between [111] and [110] and the third for points 
between, [111] and [001]. These values of |C(k)| are in very good agreement 
for the first three rings and the first two figures are in fair agreement for 
alltherings. The deviation of the third value of |C(k)| for the rings furthest 
from the ‘111) Brillouin zone is due to the attraction of the Fermi surface 
towards the (002) Brillouin zone. In addition, to this check on the useful- 
ness of the two beam approximation we can make an estimate of the 
importance of the beams neglected by the theory. For electrons travelling 
in the [110] direction there are, by symmetry, two diffracted beams of 
equal amplitude. The two beam approximation holds less well for such 
electrons. than for electrons travelling in other directions since it ignores 
one of these two diffracted beams entirely. The amplitude of the other 
diffracted beam is (using eqn. (8)) approximately + of the amplitude or + 
of the intensity of the direct beam. We may argue roughly therefore 
that the diffracted beams ignored in the two beam approximation all 
have intensities less than 3; of the intensity of the direct beam. 


§ 3. A CALCULATION OF THE RESISTIVITY OF STACKING Fauuts In TERMS 
oF ELecrron DrirrRacTion 


3.1. The Effect of Stacking Faults on the Diffracted Beams 


Seeger (1956) investigated the scattering in a crystal containing a 
fault as shown in fig. 1. He used a perturbation potential V'(r). 
V'(r)=0 in region I, 
V'(r)= —V(r)+ V(r+R) in region IT 
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where V(r) is the potential in the perfect crystal and R is the displacement 
vector associated with the fault. Since V’’(r) is either zero or has the same 
periodicity as the lattice the scattering will be from states of wave vector k 
to those of wave vector k’ =k + g where (1/27)g is a reciprocal lattice vector. 
Perturbation theory therefore demonstrates the intimate connection 
between the scattering from the fault and the normal process of Bragg 
diffraction in the crystal. However, there is no justification for its further 
use since the perturbation potential V’(r) is as large as the crystal 
potential itself. 


Reflected Transmitted 
Wave 


Schematic diagram of a stacking fault which lies on a (111) plane and produces 
a shear so that the crystal in region II is translated relative to the crystal 
in region I by a constant amount R=(a/6) [112]. The incident, reflected 
and transmitted waves used in the wave matching calculation are shown, 
but not the Tamm waves. 


We can obtain a more precise description of the effect of the fault by 
noting that if the crystal potential in region I (fig. 1) is V(r) (eqn. (4)) 
and the crystal potential in Phe II is V(r +R) then 


V(r+R)= a3. U, exp [1g .(r+R)]. 


It then follows that in region I C(k) is given by eqn. (7) and in region IT 

U. eoxp (ig .R) (9) 
H(k+g)?—k3} 4+ V/{[H(k+ 8)? —k P+ [U,P} 
Thus we see immediately that the effect of the fault is to leave the 


relative amplitudes of the direct and diffracted beams unchanged but 
to alter the relative phases of the two beams by an angle a=g.R on one 


C(k) = — 
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side of the fault. It is worth noting that the importance of this phase 
effect for the scattering of high energy (100kv) electrons in thin foils 
containing stacking faults has already been realized (Whelan and Hirsch 
1957) and applied to explain the contrast effects in the electron microscope 
image due to such specimens. The phase angle « can be shown to take 
the values 0°, +120° depending on the diffracting planes involved. If 
these planes are parallel to the fault g and R are mutually perpendicular 
so that ~«=0 and the fault has no effect. This point has been fully appre- 
ciated in the investigations made of the scattering effect of stacking faults 
on, X-rays (Paterson 1952) and high energy electrons (Whelan and Hirsch 
1957). However, it does not appear to be noted in any of the literature 
dealing with the scattering of conduction electrons by stacking faults. 
We see immediately that a fault on a (111) plane will have no effect on 
conduction electrons travelling in the [111] direction since these electrons 
will be diffracted by planes parallel to the fault. The electrons most 
strongly affected by the fault will be those travelling in the [111], [111] 
and [111] directions. In Seeger’s calculation in which a spherical Fermi 
surface was used many of the diffraction effects are averaged over angles 
so that the simple physical picture outlined above is obscured. 

Equations (7) and (9) describe the equilibrium situation between the 
direct and diffracted crystal waves at large distances on either side of the 
fault and show that is it only the phase relation between the two beams that 
is different on opposite sides of the fault. The detailed way in which the 
crystal waves change from one equilibrium situation to the other as the 
electron passes through the fault involves multiple Bragg scattering and 
interference between direct and diffracted beams from either side of the 
fault. Without, however, inquiring closely into these processes whose 
study belongs to the dynamical theory of diffraction, it is possible to come 
to two important conclusions. 

Firstly we see that the distance required to set up the new equilibrium 
will depend on how strongly the direct beam is diffracted. According to 
the dynamical theory the direct beam (when not maintained by an applied 
field) will be completely diffracted after travelling a distance known as 
the extinction distance ft) given by t)=AH/|V,| (Heidenreich 1949). The 
extinction distance is therefore a good measure of the distance travelled 
by the electron before equilibrium is restored. For conduction electrons 
in copper diffracting from (111) planes the extinction distance is of the order 
of the lattice spacing. 

The second conclusion, is that a knowledge of what the final result of the 
scattering must be is sufficient to derive an expression for it. We make 
the assumption that the direct beam is maintained constant by the field 
and thus suffers no depletion as a result of the scattering. The diffracted 
beams, however, are purely secondary effects due to the presence of the 
crystal and their amplitudes and phases relative to the direct beam on 
either side of the fault are determined by eqns. (7) and (9) which express 
the conditions for dynamical equilibrium between the direct and diffracted 
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beamsin each region. We assert, therefore, that irrespective of the detailed 
mechanism of the scattering, the nett effect is that enough additional 
electrons are scattered from the direct to the diffracted beam to set up 
equilibrium again. The additional intensity scattered into the diffracted 
beam is then 
A(k, k +g) = |O(k) exp (ig .R)—1P 
= 4| C(k)[? sin? da. Mee ee 6110) 


Section of the Fermi surface in copper as determined by Pippard. The surface 
is divided into rings near the [111] direction (see text). 


It is not permissible to regard this as a scattering probability since it 
may well be greater than unity. Values of A (k,k+g) greater than unity 
would of course be physically unreasonable in the absence of an applied 
field when there is only a definite number of electrons in the direct beam 
and any scattering diminishes that number. However, in the presence of 
the applied field the number of electrons in the direct beam is maintained 
constant by the acceleration, of new electrons by the field and thus, since 
the effect of the fault extends for an infinite distance behind it, the scattering 
will continue until eqn. (10) is satisfied. 
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The rest of the calculation consists of evaluating A(k, k+g) for various 
values of k on the Fermi surface and then. of estimating the total scattering 
produced, for several directions of the applied field, by a given density 
of stacking faults equally distributed on all four sets of (111) planes. 
Table 1 gives the value of « and hence, by eqn. (8), of |C(k)| for each ring 
on the Fermi surface near the Brillouin zone (see fig. 2). The simplest 
assumptions are made about the current distribution over the surface. 


Table 1 


Ring No. ] 2 3 4 5) 6 


€ 0-025 0-075 0-125 0-175 0-225 0-275 
0-87 0-68 0-54 0-44 0-36 0-31 
|C(k)] 0-86 0-65 0-53 0-39 0-31 0-27 
0-86 0-65 0-53 0-44 0-41 0-41 
ig 0-014 0-035 0-046 0-060 0-074 0-086 
0-011 0-011 


The first row of values of |C(k)| are used in the calculation. 


Table 2 
Field Ring 
Direction No. Y1 Y2 ¥3 Ya 
1 0-30 0-30 0-30 0-30 
[100] 
2=6 0-39 0-39 0-39 0-39 
1 0:30 0 0 0-30 
[110] 
2-6 0-71 0 0 0-71 
I 0 0-20 0-20 0-20 
[111] Lae oly fe 
2-6 0-82 0-20 0-20 0-20 


The number of current carriers from each ring is taken to be proportional 
to the area of the ring projected in the direction of the applied field. The 
velocity of the electrons assumed parallel to k must also be projected in 
the field direction. The product y of these two projection factors has 
been evaluated for various rings and field directions and is given in table 2. 
Yr» Ya Y3 and y, apply to electrons near the (111), (111), (111) and (111) 
Brillouin zones respectively. The area of each ring is given in table 1 
as a fraction, f, of }S. Sis the total area of the Fermi surface and the factor 
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3 is the average value of the projection factor y=cos?¢ for a spherical 


surface. 

The faults are assumed equally distributed on, all four slip planes and 
each fault area must be projected in the direction of propagation of the 
electrons considered. This projection factor, multipled by the appropriate 
value of 4sin?a/2=0 or 3, gives a further factor of ° for each Brillouin 
zone. ‘The fraction of the current scattered for unit area of stacking fault 
evenly distributed over all four slip planes is now given by 

x=2  fIC(k)P Diy. 
rings i 

The factor 2 takes account of the four zones on the other side of the 
Fermi surface. We find that y varies by about 10% for the three field 
directions taken (being greatest for the [100] direction) and an average 
figure is y=0-15. This fraction of the current is, to a fair approximation, 
scattered to the other side of the Fermi sphere, i.e. it is reversed in sign. 
Thus for the rate of change of current due to scattering we have 


J 
(5) = —2yJA 
ot scattering 


where A is the area of fault encountered per second. If v is the Fermi 


velocity then the average velocity, el. in any direction is $v for a spherical 
Fermi surface and, if 8 is the density of stacking faults, A = 4$fv. 


oJ 
So (5) = — yfbvJ. 
ot, scattering XP 
If the rate of displacement dk of the Fermi surface by the applied field F 


is given by 8k=eF/f the current increases according to the relation 


(5) Sit | vds .3k 
ot field 4a Fermi surface 


euSF , : : 
= for a spherical Fermi surface of area S 
7 


where p is the resistivity. 


pie oJ 
In equilibrium (5) = -(5) 
: ot scattering ot field 
9-3 4 
ES 
Taking the area of the Fermi surface S= 2-2 x 10'7em~ and y=0-15 as 
found previously we obtain 


sO 


p=1-05 x 10- Bohm-cm. 


It appears therefore that the stacking fault resistivity predicted by this 
simple model is in good agreement with the existing experimental results 
(eqns. (1), (2)). Doubtless the answer would be altered if a better approxi- 
mation than the two beam approximation were used or if account were 
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taken of the variations of electron velocity over the Fermi surface but it 
is not expected that these refinements would affect the general order of 
magnitude of the result. 


§ 4. Marcuine or Brocnh Wave Funcrions at Stackine Favs 
4.1. Reflection of Electrons at a Stacking Fault 


It is possible to obtain a useful check on, the results derived in the last 
section in terms of a more formal mathematical procedure in which the 
various crystal waves that can propagate in the perfect crystal on either 
side of the fault are matched so that the wave function and its derivative 
normal to the fault are continuous across the fault. In addition to a wave 
incident on the fault in region I (fig. 1) and a transmitted wave in region II 
it is necessary to have a reflected wave travelling back into region I. 
For this reason it is possible to refer to the reflectivity of the fault although 
the details of the way in which the reflected wave arises from the processes 
of multiple diffraction which operate at the fault are not clear. Since the 
waves must match at all points on the fault plane it follows that all the 
waves used must have the same tangential components of wave vector. 
Consequently if the incident, reflected and transmitted waves each have a 
direct beam and any number m of associated diffracted beams the require- 
ments of continuity of the wave function and its derivative will yield 
2(m +1) equations (two for each component beam) to be satisfied. Since 
it is possible to choose only the amplitudes of the direct components of 
the reflected and transmitted waves the problem can only be solved in 
the case m= 0 which corresponds to the free electron approximation when 
the fault has no effect. As Seeger (1956) has pointed out it is necessary to 
employ in the matching technique not only the waves which propagate 
in, an infinite crystal but also waves localized at the surface of the fault. 
These waves will now be discussed. 


4.2. Tamm Surface Waves 


Tamm (1932) first pointed out that it is only in an infinite crystal that 
all the components of k, the wave vector of a Bloch wave function 
exp (tk. r)U, (r) need be real. Ina finite crystal there are localized surface 
states corresponding to complex components of k normal to the surface. 
Such states have energies in the forbidden bands. Later authors (for 
references see Shockley 1939) have investigated some of their properties 
but they do not appear to have been extensively used. Attention seems 
moreover to have been restricted to the case when the Bragg diffracting 
planes are parallel to the surface in question. We shall require waves 
which are attenuated in the [111] direction and suffer diffraction from 
(111), (111) or (111) planes. 

The conventional two-beam approximation used in §2.2 leads to a 
dispersion, equation, 


(hy? + ky? + keg? — y}{(ky +94)? + (hy + 95)? + ks? —y}+ |[U,P?-=0 rat (Lay 
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where y=(2m/h?)E,—U, (Ey being the Fermi energy) and the axes 
(x, y, z) have been chosen so that 


k= (ky, ky, ks) 


and $= (91; Jo, 9). 

Equation (11) is normally regarded as a quadratic equation for y whose 
solution is given, in eqn. (5). However, if we wish to match waves of known, 
energy on the plane x=0 it can be used as a quartic equation in k,, the 
component of k normal to this plane. In such wave matching calculations 
all allowed waves must have the same tangential components of k; i.e. 
k, and k3 can also be taken as fixed. In the cases that we consider two of 
the roots of eqn. (11) are real and two complex. If the roots are written 
k,, ky’, p+iq we see that k, and k,’ represent incident and reflected waves 
respectively and the two complex roots describe surface waves decaying 
in opposite directions. Before any matching calculation can, be performed 
it is therefore necessary to know the reflected wave vector k’ = (k,', k., kg) 


Fig. 3 
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[By courtesy of the Royal Society] 
Stereographic map showing contours of equal radius vector on Pippard’s 
Fermi surface. To express the distances in the units of |g| used in the 
text the figures given should be multiplied by 0-45. 
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and the Tamm wave vectors p.. = (p + 7q, ky, ks) associated with any incident 
wave vector (k,, ky, kz) and energy H, and also to know the amplitude of 
the diffracted wave connected with each one of these waves. 


Fig. 4 


t 
O4155 0-079 O -0-079 -O-I55 


Stereographic contour map for Pippard’s Fermi surface showing as solid lines 
the contours of constant ky, (component of k parallel to [211] and as 
dotted lines the contours of constant ky (component of k parallel to 
[O11]). The two points of intersection P and @ of any two contours 
correspond to a pair of incident and reflected waves. Note how these 
points do not always lie on opposite sides of the line k,=0 (the (111) 
plane shown with alternate dots and dashes). In these cases refraction 
occurs (see text). The positions of the points P, and @, used in the 
matching are shown. The point P, cannot be shown because it is on 
the other side of the Brillouin zone. It is represented by the point 
P,' which has a wave vector directly opposite to that of the diffracted 
wave associated with P,. 


The method used to find the roots of the quartic eqn. (11) was again 
based on Pippard’s model of the Fermi surface of copper. We take a 
stacking fault on a (111) plane and consider electrons moving near the 
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[111] direction. It is convenient to take axis x, y, z along [111], [211] 
and [011] respectively. With the aid of the stereographic contour map 
(fig. 3) that Pippard (1957) has given, a map showing contours of equal 
k, and equal ks in the vicinity of the [111] direction has been constructed 
(fig. 4). Ifthe Fermi surface were spherical the contours of equal k, and 
k; would be small circles about [211] and [011] respectively but they are 
distorted by the attraction of the surface to the zone boundary. A given 
pair of contours k,=constant and k,=constant intersect in two points 


Table 3 

Point Py | Fas Ps Py Bs Pe Py 

deg 0-435 0-478 0-454 0-449 0-428 0-412 0-412 

ee 0 i) 0-079 0-079 0-079 0-079 0-155 

ky 0-226 0-268 0-164 0-194 is 0-201 0-214} 0-118 

ky’ — 0-099 0-087 0-081. 0-012 “0.099 — 0-145 — 0-042 

Pp 0-270 0-155 0-211 es 0-230] 0-282 0-298 0-295 

q 0-327 0-257 0-268 0-295 0-350 0-380 0-388 

C(k) 0-845 1-50 0-833 0-910 0-726 0-650 0-500 

C(k’) 0-403 0-795 0-607 0-490 0-403 0-382 0-371 
O(pZ) — 0-070 0-264 0-253 0-124 | —0-165 | —0-435 | —0-328 


41-902 | 40-8572] 41-222 | +1-462 | +2-122 | +2-432 | +2-472 


corresponding to the incident and reflected waves. By measuring the 
angles to these points from [111] and using fig. 3 the values of k, and ky’ 
are quickly obtained. Reference to the line k,=0 (the (111) plane) in 
fig. 4 shows that one effect of the deviation of the Fermi surface from 
spherical form is that k,4 — ky’, i.e. the reflection is no longer specular. 
Indeed for some points k, and k,’ both have the same sign so that it appears 
that the direct beam is refracted rather than reflected. It is interesting 
to note, however, that the current J, which on the average flows in the 
direction, k + |O(k)|?(k +g) is in all cases reflected (though not necessarily 
specularly reflected) at the fault. After k, and k,’ have been determined by 
this method p and q (the real and imaginary parts of the Tamm wave 
vector component normal to the fault) are obtained from the relations 
between the roots of eqn. (11) 


ky +k +2p=—29,, ! (12) 
2p(ky + ky!) + hy ky! +p? +P =gy? the? + (hy + Go)? + 2(ks? — y) 


P.M. 8 
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Values of ky, ka, kg, ky’, p and q in units of |g| for several points P,, P,...P 
on the Fermi surface near the [111] direction are given in table 3. These 
points correspond to different incident wave vectors k and are shown in 
fig. 4 together with the associated points @,, Q,...Q, which correspond 
to the reflected wave vector k’. It should be noted that all the values of 
|q| are in the range 4|g| to 3[g| so that 1/|q| the “penetration depth’ of 
the Tamm waves, is rather less than a lattice spacing. Thus 1/|q|~t) the 
extinction, distance (see § 3.1). 

Values of C(k), C(k’), C(p..) and C(p_), the last two quantities describing 
the diffracted component of the Tamm waves, are also given in table 3 
and have been calculated from eqns. (6) and (7) with the origin of r so 
chosen that U, (eqns. (4), (6), (7)) is real and negative. It should be noted 
that because of the symmetry of the problem the situation is the same in 
the vicinity of the (111), (111) and (111) zones and also that the problem 
is symmetrical in + ks. 


4.3. Matching of Waves at the Fault 


The waves to be matched are two beam solutions of the Schrédinger 
equation in region I and II of the crystal shown in fig. 1. In region I 
we have: 


An incident wave 
,(r) =exp (tk. r) + C(k) exp [i(k +g).r]. 
A reflected wave 
y'(r) = Rf{exp (vk’ . r) + C(k’) exp [i(k’ +g). r]} 
and a Tamm wave 
y(r) = Afexp (ip_.r)+C(p_) exp [?(p_+g).r]}. 


Because of the fact that the current which is 7a important physical 
quantity, is always reflected at the fault the wave y,'(r) is always considered 
to propagate in region, I. 

In region II we have: 


A transmitted wave 
by(r) = T'{exp (tk. r) + C(k) exp (ia) exp [7(k +g). r]} 
and a Tamm wave 
by (r) = Biexp (ip, .r) + O(p,) exp (ia) exp [i(p, +g). rj} 


where the phase factor « is —120°. The requirement that the wave 
function be continuous at the fault gives the two equations: 


14+R+A=B+T7, 
C(k) + RC(k’) + AC(p_) = Bexp (ia)O(p,.) + 7 exp (tx) C(k). \ 
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The requirement that the normal derivative of the wave function be con- 
tinuous at the fault yields two further equations: 


ky + Bk,’ + A(p—iq)= B(p+iq)+ Thy, 
(ky + 9x)C(k) + R(ky' + gy)C(k’) + A(p—iq+ 91)C(p_) 
= B(p + iq +9,) exp (ix)C(p,) + T(ky + 91) exp (ix)C(k). 


The solutions to these equations for the real (R,, A,, B,, 7',) and imaginary 
(Ry, Ay, By, 7's) parts of the wave amplitudes have been obtained on a 
computer and are given in full in table 4. The first set of seven values 
refer to the situation when the points P,, P,,...P, in table 3 represent the 
incident beam. The second set of seven refers to the situation when 
these points represent the reflected beam and the incident beam is 
represented by the points Q,, Q, . . . Q, (see fig. 4). 


Table 4 
Point R, R, T, 7, Ae A, By B, 
Lee — 0-262 | —0-383 | +0-218 | —0-863 | —0-102 | —0-335 | —0-336 | +-0-289 
ses —0-481 | —0-417 | —0-892 | +0-069 | —0-664 | —0-239 | +.0-264 | +1-008 
Jze —0-746 | —0-406 | —0-456 | —0-401 | —0-292 0-051 —0-161 | +0-391 
Jes —0-519 | —0-455 | —0-354 | —0-725 | —0-298 | —0-252 | —0-285 | +0:457 
oe —(-235 | —0-377 | +0-383 | —0-791 | —0-077 | —0-264 | —0-283 | -+0-188 
ae —0-149 | —0-324 | +0-620 | —0-670 | —0-001 | —0-218 | —0-226 | +0-135 | 
P, —0-359 | —0-456 | +0-407 | —0-686 | —0-075 | —0-136 | —0-161 | +0-090 
(4P —0-251 | +0-369 | +0-884 | —0-113 | —0-022 | —0-192 | —0-172 | —0-280 
Qs —0-230 | +0-208 | +0-055 | —0-900 | —0-430 | —0-377 | —0-777 | —0-118 | 
Qs | 0-654 +0-353 "40-583 —0-164 | —0-047 —0-215 0-168 —0-316 | 
Q, —0-404 | +0-336 | +0-746 | —0-238 | —0-027 | —0-218 | —0-219 | —0-332 | 
Qs —0-271 | +0-436 | +0-879 | —0-004 | +.0-023 | —0-185 | —0-157 | — 0-238 
(Ay —0-196 | +0-427 40-912 +0-031 | +0-025 | —0-168 | —0-156 | —0-186 
Q7 —0-388 | +0-493 | +0-762 | +0-236 | +0-065 | —0-114 | —0-081 | —0-146 


In view of the large number of manipulations that have been, performed. 
with quantities (the C(k)’s for instance) that cannot be accurately known. 
less significance is attached to individual figures in table 4 than to the 
general order of magnitude of the results. Values of |RP=R,?+R,* 
(the second row of values referring to the second set of results in table 4) 
and of |7|?=7',?+ 7,” are given in table 5. For some reason, not at present 
understood, the values of |7'|? seem to be almost the same irrespective of 
whether k describes the incident or reflected beam and only one set is given. 


$2 
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The table shows that there is certainly a rough correlation, between high 
values of |R|? and low values of |7'|?. However, it is difficult to interpret 
the values of |R[? directly since the velocity of the reflected electrons may 
be quite different to that of the incident electrons. It should also be born 
in mind that the waves used were not normalized. Neither of these 
difficulties apply to the interpretation of the values of |Z}? which show 
immediately that over a large area (perhaps almost all) of the Fermi 
surface the ‘ transmission coefficient’ of the fault will be of the order of 
0:75 say. Moreover we may regard the fraction of the electrons not 
transmitted (25°% say) as being scattered through a fairly large angle (40° 
which is the scattering angle for specular reflection in a (111) plane from 
the [111] incident direction, might be taken as a typical angle). This is 
to be compared with the results of the diffraction calculation made in §3 
where it was shown that a fraction y= 0-15 of the electrons were scattered 
through an angle of almost 180°. On the basis of these approximate 
considerations it appears that the resistivity predicted by the wave match- 
ing calculation would be rather less (perhaps by a factor of 4) than that 
obtained from the diffraction calculation. 


Table 5 
Point dee tes bee saa 
| RI? 0-21 0-40 0-72 0-48 0-20 0-13 0-34 
| RI? 0-20 0-10 0-56 0-28 0-26 0-22 0-40 
eats 0-79 0-80 0-37 0-65 0-77 0-83 0-64 


It should be noted (see table 4) that the Tamm waves used in the wave 
matching calculation have in most cases rather small intensities. For 
the points P,, P;, P,, P, |A|? and |B|? have values of the order of 0-1 and 
often less. For the other points P,, P;, P, higher values occur the largest 
being |B|?=1-1 for P,. The reason why these higher values should be 
associated with the cases when the direct beam is refracted at the fault 
is not understood. It would of course be of interest to compute the total 
charge density bound in surface states at the fault but this involves not 
only the electrons considered so far but also those inside the Fermi stirface 
with lower energies. It seems unlikely that the average Tamm wave 
intensity for all these electrons is greater than one tenth of the intensity 
of the unattenuated waves. Since the Tamm waves appear to have a 
“penetration depth” of about one-third of a lattice distance a tentative 
estimate would be that for every thirty atoms at the fault one electron is 
bound in the surface states. As a result of this charge density at the fault 
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the wave functions would be modified in such a way as to keep the Fermi 
level constant throughout the metal. The charge density would then be 
replaced by a scattering potential. Possibly the extra scattering produced 
would become important for wide ribbons of stacking fault but for ribbon 
widths of the order of ten atomic spacings we would expect the charge to 
behave rather like a line of vacancies. The charge density quoted above 
would be equivalent to one vacancy every three atomic distances along the 
dislocation and this would increase the scattering by only a few per cent. 


§ 5. CONCLUSIONS 


The scattering of conduction electrons by stacking faults has been, 
shown to be due to a phase change in the diffraction from those crystal 
planes which do not contain the displacement vector of the fault. By 
using a simple physical model it has been possible to obtain an estimate 
of the magnitude of the effect which agrees well with experimental data 
and has been roughly confirmed by a wave matching calculation. The 
resistivity of extended dislocations in copper is therefore mainly due to 
the stacking fault ribbon. 

The two calculations are also in agreement that the scattering from two 
faults should be independent provided they are further apart than one 
lattice spacing. This means for instance that it is justifiable to compute 
the effective area of a tetrahedron of stacking faults simply by adding up 
the areas of the faces. 

It is of great importance to know whether the resistivity of stacking 
faults depends critically on the shape of the Fermi surface and thus varies 
greatly from one metal to another. We know that the resistivity will be 
_ proportional to an average value of the intensity of the diffracted beam 
for the electrons at the Fermi level but since the calculations described 
have dealt with a specific model of the Fermi surface it is not easy to draw 
any other general conclusions. However it is possible to calculate the 
stacking fault resistivity corresponding to any particular model of the 
Fermi surface. For instance, a rough estimate based on the Fermi surface 
proposed by Heine (1957) indicates that in aluminium, where the electrons 
are almost free, the stacking fault resistivity may be only about one-tenth 
of the value in copper. A low value of stacking fault resistivity might 
also be expected in nickel where the Fermi surface is a sphere containing 
0-5 electrons per atom and does not come near any of the Brillouin zones. 
However, it should be borne in mind that the two-beam approximation 
does not provide a very good description of the conduction electrons in 
such a metal. In gold and silver the Fermi surface is probably rather 
similar to that in copper although, if the energy gap at the Brillouin zones is 
much less than the value of 6-5 ev assumed for copper, the surface may not 
make contact with the zone. Ziman (1959) has calculated that if the 
energy gap at the zone in copper were reduced to 5-6ev contact would 
just occur. For such a Fermi surface the stacking fault resistivity would 
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be about three-quarters of the value calculated for Pippard’s model of the 
Fermi surface. It thus appears that the stacking fault resistivity in these 
metals will still be considerable even if the Fermi surface does not touch 
the zone but only comes very close to it. 

The wave matching technique using Tamm waves could be greatly 
improved in, power and scope if three beams could be considered instead of 
two. At present the method is restricted to regions of the Fermi surface 
where both the incident and reflected beams are diffracted by the same set 
of planes. It would be interesting to apply the technique to simple tilt 
and twist boundaries as well as to stacking faults. Such studies would give 
useful information about the resistivity of undissociated dislocations 
(composing the boundary). 

In connection, with the resistivity from undissociated dislocations it 
should be noted that there exists a diffraction effect somewhat similar to 
that due to stacking faults. A dislocation of Burger’s vector 1/2 [110] 
produces a phase change in the diffracted beam which is zero for the (111) 
and (111) reflections and 180° for the (111) and (111) reflections. This 
phase change is not, however, abrupt as in the case of the stacking fault 
(except near the centre of the dislocation) but takes place gradually so 
that the crystal waves can adapt themselves continuously to the changing 
conditions. However, if the phase change is significant over a distance 
comparable to the extinction distance the crystal waves will not be able 
to change rapidly enough and a large scattering will result. This will be 
true for electrons passing within about an extinction distance of the disloca- 
tion so that an undissociated dislocation should behave like a stacking fault 
of width one or two extinction distances. It is of interest to note that the 
close connection between, the extinction distance and the distance from 
the dislocation within which this diffraction phase change produces _ 
strong scattering is also found in the case of high energy electrons passing 
through thin foils (Hirsch et al. 1960). For these electrons the extinction 
distances are about one hundred Angstroms and the image of the dislocation 
which is observed in the electron microscope as a result of the scattering 
is approximately of this width. 

In copper this type of scattering might be responsible for about one-fifth 
of the resistivity due to the dislocations (assuming the stacking fault ribbon 
width to be ten atomic spacings and the extinction distance to be one 
atomic spacing). In aluminium where the ribbon width is only about one 
atomic spacing but the extinction distance is perhaps ten atomic spacings 
the effect ismuch more important and using these figures together with the 
estimate given above for the stacking fault resistivity in aluminium one 
can estimate a rough figure for the resistivity due to a dislocation density of 
Ncm~ in this metal of p=5 x 10-° Nohm-cm. The effect may also be of 
importance in nickel. Thus, although detailed calculations have not been 
carried out, the diffraction scattering from dislocations is expected to make 
a contribution to the resistivity several times greater than that which would 
be predicted on the Hunter—Nabarro treatment. 


Electrical Resistivity of Stacking Faults 271 


ACKNOWLEDGMENTS 


My thanks are due to Professor N. F. Mott, F.R.S., and Dr. W. H. Taylor 
for their interest and encouragement, and to Drs. P. B. Hirsch and M. J. 
Whelan for many stimulating discussions about diffraction by stacking 
faults and dislocations. I am also grateful to Dr. J. Friedel, Dr. W. M. 
Lomer, Professor A. Seeger and Dr. J. M. Ziman for valuable comments 
and to Mr. C. Hunter for assistance with the use of the high speed computer. 
Finally I wish to thank the Master and Fellows of Trinity College for a 
Research Studentship. 


REFERENCES 


ee tee J., Ham, F.8., and Korxtmr, J. 8., 1956, Bull. Amer. phys. Soc., 

wi: 

BaAvERLE, J. E., and Kornumr, J. 8., 1957, Phys. Rev., 107, 1493. 

Broom, T., 1952, Proc. phys. Soc. Lond. B, 65, 871. 

Broom, T., and Barrett, C. S., 1953, Acta Met., 1, 305. 

CHRISTIAN, J. W., and SpREADBOROUGH, J., 1956, Phil. Mag., 1, 1069. 

CLAREBROUGH, L. M., Harcreaves, M. E., and Wsst, G. W., 1955, Proc. roy. 
Soc. A, 282, 255; 1957, Acta Met., 5, 738. 

Harrison, W. A., 1958, J. Phys. Chem. Solids, 5, 44. 

HEIDENREICH, R. D., 1949, J. appl. Phys., 20, 993. 

HEIDENREICH, R. D., and SHocktey, W., 1948, Report on a Conference on 
Strength of Solids (London: The Physical Society). 

Herne, V., 1957, Proc. roy. Soc. A, 240, 340. 

Hirscu, P. B., Howie, A., and WHELAN, M. J., 1960, Phil. Trans. A (to be 
published). 

Hunter, S. C., and Nasarro, F. R. N., 1953, Proc. roy. Soc. A, 220, 542. 

Kuemens, P. G., 1953, Aust. J. Phys., 6, 122; 1956, Canad. J. Phys., 34, 1212. 

Paterson, M. S., 1952, J. appl. Phys., 28, 805. 

Prpeparp, A. B., 1957, Phil. Trans. A, 250, 325. 

SEEGER, A., 1956, Canad. J. Phys., 34, 1219. 

Seecer, A., BERNER, R., and Wo tr, H., 1959, Z. Phys., 155, 247. 

Sreecer, A., and ScHoxrckK, G., 1953, Acta Met., 1, 519. 

Seecer, A., and STEHLE, H., 1956, Z. Phys., 146, 242. 

SHOCKLEY, W., 1939, Phys. Rev., 56, 317. 

Sricox, J., and Hrrscu, P. B., 1959, Phil. Mag., 4, 72. 

Srexie, H., and Szxcsr, A., 1956, Z. Phys., 146, 217. 

Tamm, I., 1932, Phys. Z. Sowjet., 1, 733. 

WHELAN, M. J., and Hrirscu, P. B., 1957, Phil. Mag., 2, 1121, 1303. 

Ziman, J. M., 1958, Nuovo Cim. Suppl. 7, 353; 1959, Proc. roy. Soc. A, 252, 63. 


(27200 


Infra-red Measurements of the Optical Constants of 
Liquid Silver; 


By J. N. Hopeson 
University College of North Staffordshire, Keele, Staffordshire 


[Received October 22, 1959] 


ABSTRACT 


The optical constants of liquid silver have been measured in the infra-red 
between wave numbers 4000 and 17000 cm7!. The results are compared 
with the Drude theory and lead to a value of about 1-1 for the effective 
number of free electrons per atom. The electrical conductivity of the 
surface layer, calculated from optical constants, differs by less than 10% from 
the bulk conductivity. Some results on the temperature variation of the 
optical constants are presented. The optical constants of liquid copper 
have also been measured at two wave numbers in the infra-red. 


§ 1. Iytropuction 


THE programme of infra-red measurements of the optical constants of 
liquid metals, which was begun with mercury, Hodgson (1959), has been 
continued with silver. Silver has been chosen as a typical monovalent 
metal. It has the lowest melting point (961°C) of the monovalent noble 
metals but a new furnace was needed to reach this temperature. The new 
furnace will maintain temperatures up to 1050°c but the heating coil 
fused when the temperature was raised above this value. Attempts to 
measure liquid copper (M.P. 1083°c) were only partly successful for this 
reason. 

The optical constants of solid silver have been previously measured 
over a wide spectrum range from infra-red to ultra-violet; references in 
Schulz (1957). The infra-red values of the dielectric constant, «, at wave 
numbers less than 104 cm-!, are in agreement with the Drude free electron 
formula assuming about 1-0 free electron per atom. At wave numbers 
above 104cm~! the values of ¢ deviate from the free electron curve due 
to inter-band transitions. The infra-red values of the conductivity, o, 
agree with the Drude free electron formula if the relaxation time is treated 
as an adjustable parameter. The required relaxation time is several times 
smaller than the value deduced from the static conductivity. The ratio 
varies for different silver surfaces and is probably connected with surface 
imperfections. 
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§ 2. APPARATUS AND FurRNACE DESIGN 


The apparatus and method used for measuring optical constants were 
- the same as previously reported (Hodgson 1959). The silver was melted 
in a furnace inside a water-cooled brass tank, filled with hydrogen at a 
pressure of about 2cmHg. The furnace consisted of a silica tube, 3 cm 
diameter by 7 cm long, wound with a heating coil of nichrome wire. The 
tube was supported in a steel frame and thermally insulated with alumina 
cement. The temperature inside the silica tube was measured with a 
chromel—alumel thermocouple. The melting point of silver was reached 
with just under | kw power dissipation in the heating coil. A steel crucible 
was used at first to contain the silver in the silica tube. Two series of 
measurements were made with this crucible but it gradually deformed. 
There was also a possibility of contamination of the liquid silver by the 
steel. A second crucible was made from hard carbon rod. This crucible 
was unaffected by liquid silver and has been used for most measurements. 
The water cooling did not extend to the baseplate of the brass tank and 
it was feared that the rubber seals of the electrodes might overheat. 
The rubber O-rings were therefore replaced by lead O-rings which have 
functioned satisfactorily. The two Pyrex windows on the tank were 
held against silicone rubber O-rings by atmospheric pressure. 


§ 3. MEASUREMENTS ON LiIQuID SILVER 


About 100 ¢ of silver were needed to fill the steel or carbon crucibles. 
The silver used was of purity greater than 99-98%. <A steel and a carbon 
crucible were filled with separate lots of silver. The curvature of the 
surface of liquid silver due to surface tension was greater than with 
mercury and the reflecting area was limited to the flat central part of the 
liquid surface. The surface was very sensitive to vibration and the 
supporting table had to be well insulated from the vibration of the vacuum 
pump. When first melted the silver had a cloudy surface but the cloudy 
areas soon moved to each end of the silver blob leaving a clear area in the 
middle. This clear area showed no scattered light even at angles close to 
the reflected beam. The impossibility of cleaning a liquid silver surface 
in situ was a disadvantage of this apparatus. A surface cleaning method 
similar to that used for mercury would be difficult to apply. 

Optical constants were measured at eight points in the infra-red 
spectrum. The experimental errors have been estimated by taking 
probable errors of the measured angles and calculating the corresponding 
errors in <« and c. In fig. 1, the errors lie within the plotted circles. The 
ribbon filament lamp used as radiation source in the infra-red measure- 
ments was replaced by a mercury lamp for visible and ultra-violet 
measurements. The angle of incidence (81-5°) was unsuitable for 
measurements on silver in this spectrum range. The accuracy of the 
results is therefore low, but much better for « than o. The experimental 
values of —« and o/c in fig. 2 have probable errors + 0-4 and + 0-2 x 104, 


respectively. 
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Four sets of infra-red measurements were made over a period of several 
months, two on silver in the steel crucible and two on silver in the carbon 
crucible. These sets of measurements show differences which will be 
described later. The points of fig. 1 represent a typical set of results. 
An attempt was made to measure the temperature variation of « and o 
over the limited temperature range available, about 100°c. 


Fig. 1 


S6— 53/4 +38 S04 Ore eee 


Log k 


Liquid silver (infra-red). Experimental points for liquid silver (Ag-—1); 
theoretical curves calculated from (3) and (4) with kp? = 5-26 x 10°, 
ky=1-51 x 10°, (i.e. Np=1-13, o)/¢=174 x 104) 


§$ 4. Discussion 


, The number of atoms per em$ (V ,) and the static electrical conductivity 
in €.8.U. (09) for liquid silver, have been calculated from data in Metals 
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fieference Book (Smithells 1955). At a temperature ¢°c between the 
melting point and 1100°c: 

N,=5-20 x 10?2[1—9 x 10-5(f-961)], 3. 2. (DY 
Gae@— 162 x 10°/1=orerel0-(r=961)), =. ©. (2) 
The ratio of the conductivity of the solid to the conductivity of the liquid 
at the melting point is 2-1. This agrees with the theory of Mott (1934) 


who assumed no change in the number of conduction electrons per atom 
when a metal melts. 


Sere OAS era A 


Log k 


Liquid silver (visible and ultra-violet). Experimental points for liquid silver; 
theoretical curves continued from fig. 1. 


The experimental values of ¢«, the dielectric constant, and o, the 
conductivity in e.s.u., have been compared with Drude’s free electron 


formulae : 
1—e=k,?/(k? + kp”), i eee nee! 1 (3) 


PiCe tenet bee ee (4) 


kis the wave number in cm-!; k,2=(Ne?)/(7m*c?), where N is the number 
of free electrons per cm? of charge — e and effective mass m*; ky =1/(27¢7), 
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where + is the relaxation time. The effective number of free electrons per 


atom is defined by: 
No= (75) /(Gt)- Ae 
m* m 


oo/C=$ky"/kp- EVLA OT Ay” ee ee) 
The formula (3) for « does not include a small term due to ion core 
polarization. Using the value of polarizability of the Ag* ion given by 
Van Vleck (1932), the ion core contribution to « has been estimated as 1-1. 
Values of k,2 and k, were calculated by means of (3) and (4) for each 
pair of experimental values of « and o at the same wave number, k. For 
iog k between 3-6 and 4-0, these values of ky? and k, are constant, within 
experimental error, for each set of measurements. This may be seen in 
fig. 1 where the experimental values Ag—1 are compared with curves 
calculated from formulae (3) and (4). Average values of k,? and k, for 
each set of measurements are given in the table. Ag—l and Ag—2 refer 
to liquid silver in a steel crucible; Ag—3 and Ag—4 refer to liquid silver 
in a carbon crucible. The temperatures used for Ag-l and Ag—2 are 
uncertain because of a fault in the thermocouple circuit. Ny and o,/c 
were calculated by means of (5) and (6). The reason for the variations 
in k,? and k, for different sets of measurements has not been found and 
requires further investigation. 


When k=0, then: 


Experimental Values for Liquid Silver 


ee kp NV, O,/C t 
Ag-l 5:26 1-51 1-13 174 — 
Ag-2 5-44 1-66 1:17 164 aS 
Ag-3 5-02 1:53 1-09 164 1020 
Ag—4 5:34 1-44 1:15 185 1000 
<x 109 em-? |x 103 em— SO? °C 


The infra-red values of oo/e do not differ by more than 10° from the 
electrical values given by formula (2). This indicates that for liquid silver 
the surface layer has nearly the same conductivity as the bulk metal. 
The values of Ny for liquid silver are larger than the accepted optical value 
for solid silver, V)= 1-03 (m*/m=0-97), quoted by Schulz (1957). This 
value is probably uncertain by at least +3%. The experimental points 
at wave numbers greater than 104cm-! deviate from the infra-red 
theoretical curves as shown in fig. 2. 

The temperature variation of the optical constants was measured at 
logk=3-748 and 4-020. The constants were measured initially at about 
1000°c, then at about 1090°c, and finally again at about 1000°c. The 
initial and final values agreed within experimental error and the measure- 
ments were repeated with consistent results. Values of ko’, kp, and 
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o,/¢ (infra-red) were calculated from the experimental values of ¢ and o. 
The temperature coefficients have been calculated as fractions of the 
initial values (i.e. near the melting point). 


Parameter keg? kp o,/C 
Temperature coefficient per °c = —2x 10-4 + 6:x 10-4 28 10"* 
(probable error +1 x 10-*) 


an See, . 
The coefficients of ko” and o,/¢ are of the same order as the coefficients 
of V, and o,/e (electrical) in formulae (1) and (2). 


$5. MEASUREMENTS oN Liqurp COPPER 
Pieces of sheet copper of 99-998°, purity were melted in a carbon, 
crucible. The furnace maintained a temperature sufficient to melt copper 
for only a few minutes before the heating coil fused. Measurements of the 
optical constants of liquid copper were obtained at two wave numbers. 


log k —e ale kp hae 
3-645 209 30-7 2-94 5-91 
3-748 146 22°3 3:05 5:97 

x 10? x1l0?em x109 cm-? 


The average values of k)” and ky, lead to Ny=0-84 and o,/e (infra- 
red) = 99 x 104, compared with o,/c (electrical) = 141 x 10* at the melting 
point. 


§ 6. CONCLUSION 


The infra-red optical constants of liquid silver can be explained by the 
simple Drude free electron theory. The number of free electrons per atom 
is about 1-1 and their relaxation time about 7x 10-1® sec. The value of 
the static electrical conductivity derived from the infra-red measurements 
agrees fairly well with the electrically measured value. These results are 
in contrast to those for liquid mercury, which cannot be represented by 
the free electron formulae. . 
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ABSTRACT 


Residual electrical resistivity has been used as an index of the damage 
present in tungsten irradiated in a nuclear reactor at 4°K and 77°K. A 
comparison was made between four different types of specimen prepared 
from cold-drawn wire respectively annealed before irradiation at 2000, 
1500, 800°c and not at all. A fifth type of specimen was prepared from wire 
in a different state of purity. 

Four stages of recovery have been identified, the first occurring below 
—170°c, the second between —170 and 350°c with an activation energy 
rising from 0-25 to 1-7 ev, the third between 350 and 450°c with a single 
activation energy of 1-7ev and the fourth above 450°c. Second-stage 
recovery was enhanced by the presence of cold-work, affected by impurities 
and in well-annealed samples after light irradiations it was suppressed. 

A comparison of recovery with that in other metals leads to the conclusion 
that vacancies in tungsten migrate at 400°C with an activation energy of 1-7 ev. 
The general suggestion is advanced that vacancy migration occurs in all metals 
at a temperature approximately 20% of the absolute melting point 7, °K, with 
an activation energy roughly 5 x 10~4 7’, ev. 

Recovery in tungsten above —170°c is attributed to the release of inter- 
stitials from traps associated with impurity atoms and dislocation lines. 
The free interstitial is assumed to become mobile below —170°c. 

The resistance versus dose curve for cold-worked tungsten shows acceleration 
and this is considered as evidence for the action of spikes on the unstable 
arrays of vacancies produced by deformation. 

A quantitative theory of damage and recovery developed from the above 
model predicts the observed behaviour of the five types of tungsten specimen. 


§ 1. INTRODUCTION 


In this paper experiments are presented in which neutron irradiation 
damage has been studied by observing changes in the electrical resistivity 
of tungsten. It has previously been shown that after irradiation at 
— 196°c this damage recovers by thermal annealing in two main stages at 
— 80°C and at 350°c, the activation energies associated with these stages 
being 0-5ev and 1-7ev (Kinchin and Thompson 1958). Later work 
showed that the low-temperature recovery was not observed after very 
short irradiations and this led to the suggestion that low-temperature 
recovery was due to a mobile defect being released from trapping centres 
(‘Thompson 1958 a). 
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The following experiments have been designed to determine both the 
nature of the mobile defect and of its traps: 
(i) The effect of neutron dose on recovery. 
(ui) A comparison of annealed and cold worked tungsten. 
(iii) A comparison of material from different manufacturers. 
(iv) Recovery after irradiation in liquid helium. 


§ 2. EXPERIMENTAL 
2.1. Specimens 
In table 1 the various types of specimen which have been used are listed 
with some physical properties. Specimens made of Johnson Matthey 
tungsten annealed at 1500°c were used for reference in all comparative 
experiments. Specimens made direct from manufactured wire were 
found to be in a sufficiently cold-worked state for experiment (ii). 


Table 1. Specimens 


Type A B C D E 
ates : (i) Gi). e = ie 
Experiment (i) Gay aa) (11) (11) (iii) 
Tungsten 
Manufacture Johnson-Matthey Co., London Mfg. Co. 
London 
Spec. analysis 0-005% Fe spec. pure 
Temp. of pre- 
irradiation 2000 1500 500 ae 1500 
anneal, °c 
Grain size, needles needles 
microns ay 2 1 x 25 1x25 my 
Residual 
resistivity 0-04 0-06 0-30 0-36 0-09 


micro-ohm cm 


The specimens were made from 0-008in. diameter wire which had been 
cleaned electrolytically in caustic soda solution using a nickel anode. A 
7cm length of wire had potential leads of the same wire spot welded at a 
distance of lem from either end. The resulting specimen was spot welded. 
to leads passing through a pressed glass seal which was subsequently fused 
into a lead glass envelope using a purified argon atmosphere to prevent 
oxidation. The envelope was evacuated and outgassed at 350°c after 
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which type A, B and E specimens were given a pre-irradiation anneal by 
passing a current before finally sealing off. In the case of type C speci- 
mens the pre-irradiation anneal was given to the wire before mounting. 


2.2. Resistance Measurements 


All resistance measurements were carried out in liquid helium using a 
potentiometer to compare the voltage drop across the specimen with that 
across a standard resistance carrying the same current. In order to 
eliminate drifts a regulating device (Thompson 1958b) held this current 
constant to 0-001%. The residual resistance of the specimen was measured 
to an accuracy of 0-:003°%. Because of uncertainties in the dimensions 
absolute resistivity could only be measured to a few per cent. 


2.3. Irradiations 


Irradiations were carried out in three vertical holes of the BEPO reactor 
at Harwell. These holes, TE.1, TE.3 and TE.13, penetrate the reactor 
core transverse to the fuel channels and have approximately the same 
thermal neutron flux of 10!* neutronscm~*sec-!. Using sulphur detectors 
(Martin and Stanwix 1956) alongside the resistance specimens, the flux 
of neutrons with energies in excess of 1 Mev was measured and found to be 
3x 10% in TE.1 and TE.13 but 1-2 x 101! in TE.3. However, in spite of 
these variations the resistance increase for a given time and a given type 
of specimen, was constant to within 10° in any of the holes used. This 
implies that damage is chiefly due to neutrons below 1 Mev in these holes. 
In view of this constant resistance increase all irradiation doses have been 
specified as times rather than integrated neutron fluxes. 

Short irradiations in liquid nitrogen were carried out in an open lead 
glass dewar which was lowered into the TE.13 hole whilst the reactor ran 
at full power. The 12-hour irradiation was carried out in a liquid nitrogen 
cryostat (Thompson and Jefferson-Loveday 1958) in hole TE.1.  Irradi- 
‘ations lasting 1 week were performed with the nitrogen cryostat in hole 
TH.3. After all these irradiations the specimens were removed from the 
pile im liquid nitrogen and mounted in the laboratory cryostat without 
warming up. 

An irradiation lasting 85 min was carried out in liquid helium using two 
concentric dewars, the outer one containing liquid nitrogen and the inner 
one about 200 cm of liquid helium. This assembly was lowered into the 
TE.13 hole, a carbon resistance thermometer indicating when two-thirds 
of the liquid helium had boiled away. At this stage the dewars were 
removed from the pile and topped up with liquid. After waiting an hour 
for radioactivity to decay the complete assembly was removed to the 
laboratory where the residual resistance of the specimen was measured in 
situ through leads which had been attached before irradiation. 
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2.4. Annealing After Irradiation 

After irradiation in liquid nitrogen the specimens were attached to 
measuring leads and placed in a cryostat-furnace assembly mounted above 
a liquid helium bath (Stubbs and Thompson 1958). This apparatus 
enabled resistance measurements to be made in liquid helium between 
30min annealing periods at successively higher temperatures in the range 
—196 to 200°c. For annealing at higher temperatures specimens were 
removed to a separate furnace. Temperature was in all cases measured 
and controlled to +0-1°c. 

Recovery curves were obtained by plotting (R,,—R,)/(R,;—R,). Where 
fis the residual resistance after annealing at temperature 7’ and R, and 
f, are the values before and after irradiation. 


§ 3. RESULTS 
3.1. The Effect of Neutron Dose 

In fig. 1 the recovery curves for a Johnson-Matthey specimen annealed 
before irradiation at 1500°c and irradiated at —196°c for 1, 12 and 140 
hours are compared. Before each successive irradiation current was 
passed through the specimen to heat it to 150°c and to restore its residual 
resistance to Ry. To check reproducibility two other specimens were 
irradiated for 1 hour and 140 hours, their results are also plotted in fig. 1. 
The increase of resistivity of the first specimen with irradiation is shown 
in fig. 4 (0). 

Apparently the effect of decreasing dose is to reduce the proportion of 
recovery at low temperatures. The increase of resistance with irradiation 
time shows slight deceleration. 


3.2. Comparison of Annealed and Cold-worked Tunasten 
The recovery curves for Johnson-Matthey tungsten wire in the cold- 
worked, 500°c annealed, 1500°c annealed and 2000°c annealed states are 
compared in fig. 2 after 140 hour irradiations at —196°c. In fig. 3 the 
comparison is made after a 1-hour irradiation. The increase of resistivity 
with irradiation time is shown in fig. 4 for the four types of specimen. 
The presence of cold-work evidently introduced its own recovery above 
250°c (which is in agreement with the findings of Schultz (1959)) and en- 
hanced recovery in irradiated specimens below 30°c. Although there is 
little difference between the recovery of 1500°c annealed and 2000°c 
annealed specimens those which were annealed at 500°c before irradiation 
showed rather less recovery in the range 350 to 450°c. The resistance 
of cold-worked specimens increased more slowly with irradiation time 
than that of annealed specimens and showed acceleration rather than 
deceleration. 
3.3. Comparison of Material from Two Manufacturers 


In view of the difficulty experienced in obtaining tungsten with controlled 
impurities, it was decided to compare two materials from different manu- 
facturers. In fig. 5 the recovery of tungsten supplied by Johnson-Matthey 
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Fig. 3 
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is compared with that from the Tungsten Manufacturing Co. after 
both specimens had been irradiated together in liquid nitrogen for 
lL hour. Both were annealed at 1500°c before irradiation. Whereas the 
Johnson-Matthey specimen increased by 5-4 x 10~*micro -ohm-cm the 
other increased by only 1-7 x 10-4 micro-ohm-cm in the same time. From 
the residual resistivities, given in table 1, it is clear that this second speci- 
men, was in a different state of purity from the Johnson-Matthey specimen. 
The effect of these impurities is evidently to enhance recovery after short 
irradiations, particularly above 0°c. 
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Comparison of recovery in materials of different purity but the same pre- 


irradiation annealing treatment; specimens of type B and E irradiated 
for 1 hour at —196°c. 


3.4, Recovery after Irradiation at 4°K 


After 8} min irradiation in liquid helium the resistivity of a Johnson- 
Matthey specimen annealed before irradiation at 1500°c increased by 
0-0028 micro-ohm-cm. Of this increase 30° annealed after 30min at 
— 196°c and a further 20% after 30min at 25°c. Since it is known that 
for short irradiations little recovery occurs in the range —120 to 100°c 
the recovery curve may be sketched as in fig. 6. 


§ 4. Discussion 
4.1. The Resistivity Increase 


The observed change of resistivity with irradiation time may be inter- 
preted using a model based on the following three postulates. 
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(i) Irradiation produces interstitial atoms and vacant lattice sites in 
equal numbers, the former defect being mobile at liquid nitrogen tem- 
perature. An attempt to justify this assumption will be made in § 4.2. 

(ii) It also produces ‘spike’ regions of temporary local disorder con- 
taining of the order 10° atoms. The existence of such regions has been 
demonstrated in the III-V compounds by Gonser and Okkerse (1958). 
In tungsten the recoiling atom from a neutron collision is in the energy 
range (10*ev) where it behaves like a hard sphere with a mean free path 
of a few angstrom units. It therefore enters the spike region straight 
away and we may expect that all interstitials and vacancies are left as 
debris after the spikes have subsided. Thus it is unlikely that further 
spikes will affect the vacancies already present. 


Fig. 6 
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Recovery after an irradiation of 8min at 4°k; B type specimen. 


(iii) In a perfect infinite crystal at a temperature where only inter- 
stitials are mobile mutual annihilation would prevent the accumulation of 
defects. In a real specimen however, in addition to vacancies there will 
be the following sinks for interstitials: surfaces and suitable grain boun- 
daries, dislocations and foreign atoms of favourable size. These last 
two may only be temporary sinks, or traps, from which the interstitial is 
activated by heating the specimen (Lomer and Cottrell 1955, Blewitt e¢ al. 
1957). It follows, therefore, that the number of vacancies present at any 
stage is equal to the number of interstitials which have either been trapped 
on temporary sinks or been lost to the permanent sinks. 
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Let C be the total concentration per atom of interstitials or vacancies 
produced by irradiation, 7’, be the initial concentration of the rth type of 
trap, C,, the concentration per atom of interstitials on the rth trap and C, 
the concentration lost to boundary sinks of which there is a constant con- 
centration B. B must be regarded as an effective concentration which 
would apply if the sinks were uniformly distributed through the crystal. 
Consider the situation as C increases to C+dC. The extra dC interstitials 
have the choice of (7',— C,,) unfilled traps of type r, of any other unfilled trap, 
of finding a boundary sink, or of annihilating the vacancies present at a 
concentration (>C,+C,). It follows from simple proportion that the 
increase in the concentration trapped on the rth trap is: 


ae (7,—C,)dC : 
B+ S(f,—C,) + 20, +G, ’ 

dC, _ _(2,~@,) __ (1) 

iG (BL GEST,) 0 ee 
Similarly 

Ci (2) 

i ~ BXOAST ea he ee e 

7) 1/2 

rence Oy=(B+37,)1| P+ ee | -1} ees 


Solving (1) we obtain for C,, 


C,=1,| 1-exp(- 58) |. eae ee ee 


Now if the resistivity change due to interstitials being trapped on the 
rth trap is a,C,, and that due to a concentration V of vacancies is bV, then 
we may write for the resistivity increase : 


Ap= Ya,C.+ >bC.+b0, 


or Ap=b{| 1~exp(— 5) ]2,(1+ 4) +es}. ees 


In this expression the first set of terms, arising from the trapped inter- 
stitials and their associated vacancies, saturate with the same period as C 
increases. The term bC, which is due to the vacancies left behind by inter- 
stitials lost to boundaries, shows slight deceleration without saturation. 

Let us apply this analysis to predict the resistivity increase of an annealed 
sample of tungsten. First the relation, between C and irradiation time 
may be established as: 
NO oo ins oe oe 
Where n is the number of interstitial-vacancy pairs per neutron hit, o is 
the total fast neutron cross section (~10-*em?2), ¢ is the fast neutron 
flux (~3 x 10!'neutron cm sec~), and 7 is the irradiation time. Taking 
m as 30 the value of C is of the order 10-4 for a 140-hour irradiation. 
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The concentration of traps, ST, may be estimated as follows; the 
residual resistivity of the tungsten suggests an impurity concentration of 
10~, assuming the reasonable value of 500 microhmecm per unit con- 
centration of impurities. This concentration is not in contradiction with 
the spectrographic analysis. Those traps which are due to impurities 
will have this order of concentration. Dislocations might also provide 
traps, a density of 101° lines cm~ producing a concentration of order 10-4 
if ten trapping sites exist in each atomic plane, pierced by a dislocation. 
It would appear then that 10-4 is a reasonable value for ‘7’. 


Fig. 7 
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Resistance increase as a function of irradiation time calculated from theory for 
A, B, C and D type specimens with the sink concentrations given in 
table 2. 


The value of B, the effective concentration of permanent sinks, will - 
depend primarily on grain size. An order of magnitude may be obtained 
as the ratio of lattice sites at the grain surface to the total of sites in the grain, 


i.e. Br 6/dA8 


where d is the linear grain size and A the density of atoms. For 10 micron 
grains B is of the order of 10~*. 

In previous work (Kinchin and Thompson 1958) a value of 2 10° 
microhm cm was deduced for the resistivity increase in molybdenum per 
unit concentration of interstitials and vacancies. For the purposes of the 
present order of magnitude calculation let us take a,=b = 10° microhm cm 
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for tungsten which behaves similarly to molybdenum. We may then 
sketch Ap as a function of C for specimen types A, B and C as in fig. 7 using 
the estimates of table 2. 


Table 2. Estimated Sink Concentrations ( x 10+) for Four 
Types of Specimen 


Type B pase? T, Vo Temp. of 
annealing °C 
A 0-05 0-20 0-05 0 2000 
B 0-50 0-30 0-15 0 1500 
C 5-00 1-00 = 0 500 
D 5-00 1-00 — 2-00 None 


A comparison of the curves in fig. 7 with the results of fig. 4 shows that 
the theory explains the qualitative effects and gives the correct order of 
magnitude of the resistivity change. Expressed in physical terms, an 
increase in grain size favours both the trapping process and annihilation of 
vacancies. Consequently there is a greater tendency to saturate and a 
greater initial rate of increase. 

In the publication referred to above, the change in resistivity of tungsten 
irradiated at 30°C was reported to show an initial deceleration up to 5 
weeks irradiation followed by steady rise. Here again there is agreement 
with the behaviour of the two terms in eqn. (5) the first of which might be 
expected to saturate after about five weeks’ irradiation. 

Let us now consider the resistivity increase in cold-worked tungsten. 
The acceleration with irradiation time suggests that there is a process 
occurring simultaneously with the accumulation of vacancies and trapped 
interstitials, which causes a decrease in resistivity. This process might 
be the action of spikes on unstable defect configurations produced by cold- 
work. for instance it has been suggested that vacancies are produced in 
strings during cold-work (Cottrell 1957) and the effect of a spike on such 
an array might well be to cause aggregation with an associated decrease in 
resistivity. 

Suppose that such a process occurs and that one spike containing S 
atoms is produced per neutron hit. From our previous definition there 
will be n interstitials and vacancies associated with each spike. The con- 
centration of spikes is thus C/n and the concentration of lattice sites affected 
by spikes is SC/n. The effect of increasing C by dC will be to increase the 
concentration of spike sites by SdC/n. If V’ is the concentration of 
vacancies due to cold work then a concentration V’ of these extra spike 
sites will be strings of vacancies and the concentration of vacancies aggre- 
gated will be: 


aa dc. 
n 
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The concentration of cold-work vacancies will be further reduced by 
annihilation with radiation-produced interstitials. We shall neglect the 
effect of this annihilation on the aggregation process. By an argument 
analogous to that used above the following equations may be shown to 
govern the concentrations V, V’, C,, and C,: 


dV _ >(7,-—C,)+ B+ V' 
dC B+>(7,-C,)4+V4+V"’ 


dV’ SV' V’ 5 
a - s eee eee (7) 


dG ti. Baa (lO eye’ 

dO T,—O, 
Ce 2S (T= Cae 2 
aC; B 


dC B+ >(T,-C,)+V4V"’ 

with the boundary conditions 
C=03/—C,=Cr=0 

and | a rs 
A numerical solution of these equations has been obtained for the case of 
cold-worked tungsten using the values of Vy’, Band $7’, in table 2, taking 
S as 10° and n as 30. The value of V,’ was deduced from the residual 
resistivity difference of type C and D specimens. From this solution Ap 
has been calculated as a function of C and values are shown in fig. 7.. The 
predicted acceleration with C is in agreement with the observations 
recorded in fig. 6. 


4.2. Recovery 


Blewitt (private communication) has shown that tungsten irradiated at 
20°K for 1 week in the Oak Ridge pile recovers about 50° of its resistivity 
increase after annealing at liquid nitrogen temperature. Since this is in 
substantial agreement with the tungsten irradiated at low temperatures 
for 8 min it appears that the effect of dose on this part of the recovery 
curve is small. An approximate recovery curve between 4°K and 500°C 
for a 140-hour irradiation may then be sketched as in fig. 8. Recovery 
below — 170° has been designated stage I, that between — 170 and 350°C 
as stage II, that between 350 and 450°c as stage III and above 450°c 
as stage IV. 

Previous work (Kinchin and Thompson 1958) has shown that 
activation energies in stage II rise from 0-25 to 1-7ev and that stage III 
proceeds at 1-7ev. 

Recovery curves of this type have been observed in other metals, the 
difference being that in some cases, notably gold, stage I is absent (Cooper 
et al. 1954) and that in the case of electron irradiation stage II is greatly 
reduced (Corbett et al. 1957). However, in all reported cases of irradiated, 
quenched or cold-worked metals stages III occurs at an absolute tem- 
perature roughly 20% that of the melting point. Table 3 lists some 
typical activation energies and temperatures for stage III in Ag, Al, Au, 
Cu, Mo, Nb, Ni, Pt. and W. 


290 M. W. Thompson on the 


It will be noticed that the ratio of activation energy to absolute melting 
temperature varies little with the type of damage in a given metal or from 
element to element. It has long been established that such a relation 
exists between self diffusion activation energies (Kidson and Ross 1957) 
where vacancy formation and migration is the common mechanism. This 
suggests that the same stage III recovery mechanism operates after 
irradiation, cold work or quenching and that this mechanism is the same 
for all the above metals. 


Fig. 8 
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Since stage III recovery after quenching is generally attributed to the 
migration of vacancies (Koehler ef al. 1957) it is assumed that this is also 
the case in irradiated and cold-worked material. In particular, let us 
attribute stage III in tungsten to vacancy migration. 

In order to explain the low temperature recovery it is necessary to 
postulate that interstitial migration occurs either in stage I or at some 
even lower temperature. For the purposes of this discussion let us 
assign stage I to interstitial migration. 

Specimen types C and D have the same grain size and impurities and 
should have roughly the same dislocation density. Type D, however, 
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Table 3 
: E Gi 
Metal | 7',,°K | Damage| E ev hse lees 10 7 Reference 
™ m 
fe | 1234 | ow. 0-65 | 230 | 5:3 |0-17 | Manintveld 1952 
Trrad. 0-55 190 5-9 | 0-20 | McReynolds et.al. 
1955 
Al Quench | 0-44 300 4-7 0-32 | Bradshaw and 
ORE 933 Pearson 1957 b 
Quench | 0-52 270 5:6 0-30 | de Sorbo and 
Turnbull 1959 
C.W. 0-69 290 5:2 0-22 | Manintveld 1952 
Au Quench | 0-82 310 6:1 0-23 | Baurle and 
f.c.c. 1336 Koehler 1957 
Quench | 0-68 370 5-1 0-28 | Bradshaw and 
Pearson 1957 a 
| Cu 1356 Irrad. 0-68 250 5-0 0-19 | Overhauser 1955 
f.c.c a LCLWe 0-67 230 5-0 |0-17 | Eggleston 1952 
M Irrad. 1:3 420 4-5 0-15 | Kinchin and 
b 4 2893 Thompson 1958 
ee C.W. 1-26 430 4-4 | 0-15 | Martin 1957 
NP. | 2220 |Irrad. | 1:3 420 | 5-9 |0-19 | Makin 1959 
.C.C 
Sn LvZ84.0.W; 1-0 500 | 5-8 | 0-28 | Nicholas 1955 
Irrad. ]-2 420 5-9 |0-21 | Dugdale 1952 
2047 | Quench 1-1 570 5-4 |0-29 | Bradshaw and 
f.c.c. Pearson 1956 
Trrad. 1-7 670 4:7 0-18 | Kinchin and 
W 3643 Thompson 1958 
b.c.c. C.W. 1-7 620 4-7 | 0-17 | Schultz 1959 


has an initial concentration of vacancies formed during deformation. 
The recovery curves of fig. 2 show that the effect of these vacancies is to 
enhance stage II recovery. This is evidence that stage IT is in some way 
associated with interstitial migration resulting in vacancy annihilation. 
It was pointed out in an earlier publication (Thompson 1958 a) that the 
absence of stage II in lightly irradiated tungsten may be explained by 
assuming that this stage is due to the thermal activation of interstitials 
from traps in the lattice. After a one hour irradiation C'~ 10 whereas 
>7T,~ 10-4. Under these conditions interstitials being activated from 
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the shallowest traps will become attached to deeper traps rather than. 
annihilating the small concentration of vacancies present. Provided that 
the resistivity change associated with those interstitials which reach 
boundaries is balanced by that due to interstitials changing traps there 
will be no net recovery observed. In the case of type E specimens, which 
have a different set of impurities, the presence of stage II recovery after 
light irradiations might indicate that the decrease in resistivity due to 
migration to boundaries is not balanced by an increase due to changing 
traps. If the trapping mechanism of stage II recovery is accepted, then 
this difference in behaviour between types B and E is clear evidence that 
impurity atoms are capable of trapping interstitials. A similar conclusion 
has been reached by Blewitt et al. (1957) and Martin (1959, private com- 
munication) during studies of neutron-irradiated doped copper. The 
role of dislocations in recovery is not obvious from the experimental 
results, owing to simultaneous changes in grain size, but will be considered 
in the analysis which follows. 

Using the three postulates of § 4-1 let us develop a quantitative theory of 
stage II recovery for comparison with the experimental observations. In 
order to simplify the analysis the following substitutions will be made: 


Pee: 

tar Rg 
ST 

=| z 

t = BR 


H=1-—exp(—2). 


It follows from eqns. (3) and (4) that after irradiation the following 
concentrations of sinks for interstitials will be present : 


Vacancies : B(x+(t—1)£), 
Unfilled traps : B(t—1)1—£), 
Boundary sinks: B; 


.. total concentration of sinks = B(t#+ 2). 


N ow suppose that the temperature of the specimen is raised sufficiently 
to liberate the interstitials from the first trap (7',) and that these constitute 
a small proportion of the total number of trapped interstitials, i.e. 


2 ee TS 
The total concentration of released interstitials is found from eqn. (4) as 
ab(t—1)E. 


Then the concentration of these which annihilate vacancies is 


<Be—1ya| 0") 
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the concentration going to the rth unfilled trap (r4 1) is 


E)(t,—1) T 
1)H fol 
~ya| SA t+x | pap 


and the concentration lost to boundaries is 


wB(t—1)E lal . 


Now (a, +5) is the resistivity change due to unit concentration of released 
interstitials annihilating vacancies. a, is that due to unit concentration 
migrating to boundaries and (a,—a,) is that due to unit concentration 
becoming trapped on the rth type of trap. The net recovery of resistivity 
is then 


aBb(t—1)E 
+ 2(t— Y)(@.—4,)(1—B) +44}. 

The total change in resistivity due to irradiation is given by eqn. (5). 
We may therefore write the proportion of recovery due to interstitials 
being released from trap 1 as 
(Ap); = 

Ap 


(Ap), = — 


_ off 1)E [1 +a,/b](e+ (¢—1)£) + D(t,—1) (4 —a, |b) (1— B) + 4,/6 
t—x EY (t,—1)[1+a,/b]+x 
(7) 
From experiment it is known that for A and B type specimens recovery is 
zero for small doses, i.e. 
tp}: 
Ap 
Putting this condition into eqn. (7) we conclude that 


>0 as x->0. 


d(¢—1)(a4,-4,) = —a,, in A and B type specimens. . . . (8) 
In physical terms, the resistivity change due to released interstitials 
migrating to boundaries is compensated by that due to changing traps. 
Under these special circumstances equation (7) becomes: 

Ap) — a(t—1)B [1+ ay/b)(e+ (¢—1)B) + (a/b) ae 

Ap f+@  S@,-D[ltajole+e = | 
In fig. 9 observed values of (Ap)1/Ap for the recovery occurring between 
—120 and —40°c are plotted as a function of irradiation time for A, B, C 
and D type specimens. The expression (7a) has been made to agree with 
the type A and B experimental points at 140 hours by choosing suitable 
values for 7',. It will be seen that the remaining points then fall on the 
theoretical curves. This indicates that the theory predicts at least the 
correct functional dependence of recovery on irradiation time. The 
selected values of 7’, are shown in table 2 and it will be noticed that there 
is a factor of 3 between the values for specimen types A and B. Because 
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the coefficients a, are unknown for types C' and D and because there is no 
simplifying relation like (8), it is not possible to tabulate the expression (7) 
for these other specimens. 

Let us now consider the evidence for traps being associated with dis- 
location lines. The theoretical curves of resistance increase were obtained 
by assuming a decrease in }7', with increasing pre-irradiation annealing 
temperature. The temperatures involved strongly suggest that part of 
the trap concentration follows the dislocation density with the implication 
that dislocations provide trapping sites. The further fact that 7’, varies 
with pre-irradiation annealing temperature in a similar way, strengthens 
this suggestion. 


Fig. 9 


7 hours 


Recovery between — 120 and —40°c as a function of irradiation time: the solid 
lines represent theoretical curves which have been fitted to the observa- 
tions at s=140 using the sink concentrations given in table 2. 


4.3. Summary of Conclusions 


A quantitative theory which is consistent with the experimental facts 
may be derived from the following postulates : 
1. Interstitials and vacancies are produced by irradiation in equal 
numbers. 
2. Interstitials are mobile at 77°K, stage I recovery between 4°k and 
77°K being tentatively assigned to interstitials first becoming mobile. 
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3. Vacancies become mobile between 350°c and 450°c and this mech- 
anism is responsible for stage ITI recovery. 
4. Three types of sink compete for interstitials : 
(a) Vacancies. 
(6) Grain boundaries. 
(c) Traps. 
Thermal activation of interstitials from traps is responsible for 
stage II recovery (— 196 to 350°c). These traps are associated both 
with impurities and dislocations, the concentration depending on 
pre-irradiation annealing treatment. 
Irradiation produces spike regions of temporary disorder containing 
10° atoms. 
6. Cold-work produces vacancies in strings and the action of spikes on 
these is to cause aggregation. 


OU 
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Last of Symbols 
R,=resistance of specimen before irradiation. 
Rk, =resistance of specimen after irradiation. 
R,=resistance of specimen after annealing for 30 min at T°c. 

7 = Irradiation time. 

n=the number of interstitials and vacancies per neutron hit. 

S =the number of atoms affected by a spike. 

C =total concentration of interstitials and vacancies produced. 
C,,=total concentration of interstitials on the rth trap. 
O,,=total concentration of interstitials lost to boundaries. 

V =total concentration of vacancies present due to irradiation. 
V’ =total concentration of vacancies present due to cold work 

(initial value V,’). 
T',=total concentration of traps of the rth type. 

B=total concentration of boundary sinks. 

a, =resistivity due to unit concentration of interstitials on the rth trap. 

b=resistivity due to unit concentration of vacancies. 


«=C,/B. 
Lp=1+77/B. 
t=1+>T,/B. 
ge UA ie 


E=1-exp (-2). 
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CORRESPONDENCE 


A Possible Example of the Muon Radiative Decay 
with Internal Conversion of the y-ray 


By VALERIE MAYES 
H. H. Wills Physical Laboratory, University of Bristol 


[Received March 30, 1960] 


Ty making an area scan for 7—y-e, in a stack which had been, exposed to 
the cosmic rays at an altitude of approximately 13 kilometres, an event 
was found (PI. 33) in which three relativistic electrons escaped from the 
decay of a p-meson at rest. A similar event has been observed by 
I. I. Gurevich et al. (1959). The nature of the «.-meson was established 
from its origin and from its range, 613. (The average range for j.-mesons 
from 7—p decay in G5 emulsion is 595.) The pair and positron tracks 
dipped steeply down into the emulsion at angles to the emulsion plane of 
approximately 70° and 45° respectively. It was not possible to follow the 
electron tracks in more than one emulsion because the stack consisted of 
alternate layers of 600, emulsion and 1-6mm tungsten. The recorded 
track lengths were Le, = 188 yw, Le,=186 4, Le; = 242. The short length 
of the tracks and their steepness made energy determinations difficult. 
Apparent values for the energy of the pair and of the positron were 
obtained from scattering measurements, these gave 


E pair = 40 + 16 Mev, 
E positron=18+ 5Mev. 


From the opening angle of the pair an estimate of E pair of approximately 
25 Mev was obtained. The angles between the tracks are 


peye,~5°, eye,~ 60°, pe,e,~ 58°. 


Since the electron tracks are not coplanar and the directions of the pair 
and of the positron are not collinear, the event cannot be identified with 
either p> 3e or p>e+y with subsequent conversion of the y-ray. The 
event can be identified, however, with either (1) wt +>et+(et+e-)+v+0 
or (2) the coincidence of a converting y-ray with the decay point of a 
y-meson. The probability of radiative decay has been assessed by 
N. Tzoar and A. Klein (1958) as 10-° per x-meson and hence that of (1) 
is of the order of 10-7 per z.-meson. The probability of finding a back- 
ground electron pair within 10 cubic microns of the + decay was determined 
from the density of 7 in the emulsion and was of the order of 10~° per 
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y-meson, i.e. several orders of magnitudes smaller than the theoretical 
probability of radiative decay with internal conversion. In order to check 
that this did not underestimate the probability of coimcidence, an area 
sean for electron pairs was made. This showed that the density of pairs 
was not greater than the density of 7+ (400 per cm’) in the vicinity of 
this event. 

The probability of a positron producing an electron pair, either by 
conversion of bremsstrahlung or by trident production (i.e. direct pair 
production in the field of a nucleus), in the first few microns of its path 
after normal p decay, is of the same order of magnitude as the probability 
of radiative decay with internal conversion. In this event, however, the 
positron and electron pair were separated by approximately 60° and in 
the above two processes, at this energy, the electrons would be collimated 
within a few degrees. In the event observed by Gurevich et al. the pair 
and positron were indistinguishable, all three tracks lying within approxi- 
mately 10° of one another. 
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Direct Electron Microscopy of Thin Foils of Internally 
Oxidized Dilute Copper Alloys 


By M. F. Asnpy and G. C. Surra 
Department of Metallurgy, University of Cambridge 


[Received February 26, 1960] 


SEVERAL workers (Bollman 1956, Hirsch ef al. 1956, Nicholson et al. 
1957, Nicholson and Nutting 1958, Whelan ef al. 1957) have developed 
techniques for examining thin metal foils by direct electron transmission. 
This paper describes the application of these techniques to the examination 
of an internally oxidised copper 0-25% aluminium alloy. The authors 
have obtained electron micrographs which show the internal oxide particles 
in the copper matrix and the interaction of dislocations with these 
particles, 
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This alloy, rolled to a strip 200 microns thick, was internally oxidized 
by placing it in cuprous oxide and heating for 20 min at 950°c, producing 
a dispersion of alumina in a pure copper matrix. Electron diffraction of 
electropolished specimens, using a reflection technique, showed the 
dispersion to be of y high temperature alumina. This analysis was 
confirmed by preparing carbon extraction replicas from the specimen 
surface, and photographing diffraction patterns from them, using a 
selected area electron diffraction technique. 

Satisfactory foils of this internally oxidised strip were obtained by 
using the Bollman (1956) technique. The bath used was a methyl 
alcohol: concentrated nitric acid mixture in the proportion 2:1, at a 
potential of 9 volts and a temperature of —20°c. Figure 1 (PI. 34) is 
an electron micrograph of a foil of an unstrained specimen produced in 
this way. Many of the alumina particles, of diameter about 3004, are 
clearly defined as triangular plates or tetrahedra, and there is evidence 
of precipitation along preferred planes in the copper matrix. Pronounced 
contrast effects can be seen around the particles in the twin, where the 
lattice orientation is favourable for diffraction contrast. Two explana- 
tions in terms of strain fields are possible. Firstly, Nicholson and 
Nutting (1958) have shown that coherency strains can give rise to this 
type of strain pattern, so the alumina particles may be coherent with the 
copper matrix ; attempts to fit the alumina and copper lattices to allow 
coherency have so far proved unsatisfactory. Secondly, while the 
material is cooling after the internal oxidation process, the considerable 
difference in coefficient of thermal expansion of copper and alumina must 
produce strain fields around the particles. This second explanation seems 
the more probable one. The dark rings surrounding some of the particles 
in fig. 1 are due to carbon contamination of the foil in the microscope ; 
carbon builds up round particles which project from the surface of the 
foil. 

Other foils, also internally oxidized, were strained in tension by amounts 
varying from 0-5 to 6-0%. From among these foils, two, strained 4%, 
have been selected for illustration (figs. 2 and 3, Pl. 34). 

Dislocations are visible in fig. 2 bowed out between points where 
alumina particles prevent their movement. At first sight this pattern 
of dislocations seems consistent with Orowan’s theory of hardening of ~ 
metal by a dispersed phase (Orowan 1947). The foil is of the order of 
10004 thick, so that any one particle will not intersect all the slip planes 
in the plane of the foil. However, by measuring the distance apart 
(2R) of particles between which a single dislocation is bowing out, an 
estimate of the yield stress (c) of the material can be calculated, using the 
relationship: 

aGb 


c= — 


R 


(where « is a constant approximately equal to 3; G is the shear modulus 


300 Correspondence 


of the matrix; and b is the Burgers vector of the dislocation.) The 
mean calculated yield stress was 1-4 x 10*p.s.i. The measured yield 
stress was 1-58 x 104 p.s.i. 

Heavily jogged dislocations and triangular dislocation loops attached 
to particles can be seen in fig. 3. These loops have not collapsed onto 
the particles, but bow out from them. A stable configuration such as 
that seen in fig. 3 can be explained by supposing the loops to be prismatic, 
i.e. that all parts of the loop do not lie in the same slip plane. The 
authors have used a modification of Orowan’s theory first proposed by 
Hirsch (1957) to explain the formation of these loops. This is illustrated 
schematically in fig. 4. 


Fig. 4 


In 1 a pure edge dislocation, moving in the direction of its Burgers 
vector (b), has been held up by a set of particles intersecting the slip 
plane. As the applied stress rises, the screw components of the dis- 
location, formed when the dislocation took up the shape shown in 1, 
Cross slip, as shown in 2. These screw components are of opposite 
sign on each side of a single particle ; they attract each other, as in 3 ; 
and annihilate each other, as in 4. The process leaves a prismatic loop 
round each particle and corresponding prismatic jogs in the dislocation 
line. The loops can move on the surface of a cylinder whose axis is 
parallel to the Burgers vector, but they cannot collapse completely. 
This diagram shows the simplest form of this type of interaction. Any 
dislocation with an edge component should behave in a similar way. 
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A prediction from this mechanism is that the direction of elongation of 
the loops projected on to a slip plane will be parallel to a slip direction 
in the matrix metal. This prediction could be checked by electron 
diffraction experiments : these have not yet been tried. Inspection of 
fig. 3 shows that the loops appear to have a common direction of 
elongation. 


ACKNOWLEDGMENT 
The authors are indebted to Dr. P. B. Hirsch and Dr. A. Kelly for many 
helpful discussions and to Mr. R. 8. M. Revell for his help in obtaining the 
electron micrographs. All electron micrographs were taken on a 
Metropolitan-Vickers E.M.6. 


EFERENCES 

Botumann, W., 1956, Phys. Rev., 103, 1588. 

Hirsca, P. B., 1957, J. Inst. Met., 86, 13. 

Hrrescu, P. B., Hornz, R. W., and WueEtan, M. J., 1956, Phil. Mag., 1, 677. 

NicHotson, R. B., THomas, G., and Nurtine, J., 1957, Brit. J. appl. Phys., 
9, -25. 

NicHotson, R. B., and Nutrine, J., 1958, Phil. Mag., 3, 531. 

Ornowan, E., 1947, Discussion, Symposium on Internal Stresses (London : 
Inst. Metals). 

Tuomas, G., Nutrine, J., and Hirscu, P. B., 1957, J. Inst. Met., 86, 7. 

WHELAN, M. J., Hiescu, P. B., Horne, R. W., and Botitmann, W., 1957, 
Proc. roy. Soc., 240, 525. 


REVIEWS OF BOOKS 


Manual of Mathematical Physics. By Paun I. Ricuarps. (London : 
Pergamon, 1960.) [Pp. 486.] £5 10s. Od. 

Mosr dictionaries of mathematics or physics are useless: the topics are too 
intimately linked in logical chains; a single word, or phrase, cannot be pulled 
out and discussed in isolation; it carries with it strands of meaning connected 
with every other word in the list. This Manual is something different—it is a 
cryptic summary of the most important subjects in Mathematical Physics, more 
or less logically derived. It might be, as it were, the notes one prepares before 
an examination, reminding one of all the formulae and their relations, without 
going into all the rigorous details. The idea is a good one, and it seems to have 
been undertaken very conscientiously. Thus, the general pattern is sound, the 
style is extraordinarily easy for such a condensed work, and all is well indexed 
and complete with cross references. It is impossible to do more than check 
one’s own field; it seemed sensible though I must admit that there were a few 
typographical errors in the formulae. The level is that of the well-educated 
theoretical physicist who is not a specialist in the particular field in which he is 
seeking guidance. As a work of reference it should be very handy, in a library 
or in the private collection of the experienced research worker. J.M. Z. 


Fluid Mechanics. By L. D. Lanpat and E. M. Lirsurrz. (Pergamon Press, 
1960.) [Pp. xii+536.] £5 5s. Od. 


Tus book provides a concise and, on the whole, clear account of the funda- 
mental physical phenomena associated with the mechanics of moving fluids. 
The structure of fluids is not considered, except in so far as it determines the 
macroscopic variables of state and the various transport coefficients of mass, 
momentum and energy. It is written from the standpoint of the theoretical 
physicist, and little mention is made of many aspects of fluid mechanics related 
to phenomena of geophysical and engineering interest. Not much space is 
devoted to the detailed solution of particular problems or the presentation of 
empirical results, except where these illuminate the topics. On the other hand, 
there are accounts of subjects, such as thermal diffusion, geometrical acoustics, 
relativistic fluid mechanics, and the dynamics of superfluids, which are not 
usually included in treatises on fluid mechanics. 

The book starts with the motion of an ideal, but not necessarily uniform, 
fluid, and this is followed by accounts of the motion of a viscous fluid, turbulence, 
and boundary layers. These first four chapters are not altogether satisfactory, 
although the good features outweigh the bad. They make very interesting 
reading for the fluid dynamicist and contain some important new material, in 
particular a proof that a laminar boundary layer can separate only at a point 
where the wall stress is zero; but the arrangement of the material is somewhat 
arbitrary and occasionally the ideas are not presented clearly. There are next 
three brief chapters on some flows in which temperature gradients, diffusion, and 
surface effects, are important. These are followed by seven chapters on the 
motion of a compressible fluid or gas dynamics, which constitute a first class and 
stimulating account of the subject as it was a few years ago. The book con- 
cludes with chapters on relativistic fluid dynamics, the dynamics of superfluids, 
and fluctuations in fluid dynamics. 

This book, although not without some flaws, is a stimulating and welcome 
addition to the growing list of treatises on fluid mechanics. One serious fault, 
however, is that it is useless as a bibliography, since only the author and date 
are given for most of the original papers referred to. PoGsn: 
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Quantum Particle Dynamics. By J. McConnEtu. (North Holland Publishing 

Company, 1960.) [Pp. 266.] 40s. 

Tus second revised edition, appearing only two years after the first edition, 
is a readable and well presented account of quantum mechanics, quantum 
field theory and elementary particles. In a book of such short length certain 
topics have to be left out. This is especially noticeable in the last two chapters 
on nuclear physics and elementary particles. This may reduce the value of the 
book for those who wish to understand the modern theory of elementary 
particles. But it is still a valuable guide to the earlier stages in the development 
of this theory. 

The first half of the book is a detailed discussion of quantum mechanics. 
Chapter headings show the ground covered: Special Theory of Relativity, 
Introduction to Quantum Theory, Schrédinger’s Wave Mechanics, Heisenberg’s 
Matrix Mechanics, Angular Momentum Systems, of Particles, The Dirac 
Equation. Then a careful discussion of quantisation of the electromagnetic 
field is given, and applied to the photoelectric and Compton effects, Brems- 
strahling and Electron Pair Production. Finally in the last two chapters a 
brief consideration of nuclear forces and an even briefer one of elementary 
particles is given. The manner in which the topics are discussed makes this 
book very suitable for undergraduate level. The last two brief chapters should 
certainly whet the appetite of the interested reader for more detailed discussions 
of recent work in elementary particle physics. J: Gale 


Taylor-Memorial Manual of Advanced Undergraduate Experiments in Physics. 
American Association of Physics Teachers (Editor-in-Chief: THomas B. 
Brown.) (Reading, Mass.: Addison-Wesley, 1959.) [Pp. 550 with 241 
line diagrams.] 72s. 


THIS manual, dedicated to one of the foremost teachers of experimental physics 
in the United States, is designed primarily as an aid to those responsible for 
arranging advanced. practical courses in universities and colleges. It is there- 
fore much more a compendium of principles for the teacher than a working 
text-book for the individual student. The latter type of book is regrettably 
lacking, at least in Great Britain, and the former non-existent in the reviewer's 
experience untilnow. The nearest equivalent in the past would be the series of 
books on experimental physics by Searle, which had the advantages and dis- 
advantages of their single authorship. The new Manual is a cooperative effort, 
in which contributions from 191 individual physicists have been co-ordinated and. 
boiled down by a board of eight editors. The result is a mine of information, 
which any University physics department will do well to quarry, and not least 
those (not a few) brought up on and still anchored to * Searle’. Whilst it is true 
that the basic philosophy of experimental physics can be distilled out of almost 
any experiment, it is equally true that the individual student’s interest is likely _ 
to be more readily aroused by examples related to a modern context. Also a 
considerable number of new practical techniques have come out of modern 
physics and ought now to find a place in an undergraduate course. It will be 
helpful to many, for instance, to look through the section on non-linear 
oscillators, which includes a simple illustration of a relaxation oscillation, and 
the next section in which servomechanisms are discussed. 

All the main branches of physics are covered, the chapters on. electricity and 
magnetism and on nuclear physics being particularly extensive. The former 
includes microwave optics and the latter a section on health hazards. Inevit- 
ably in such an encyclopaedia the level and the detail of treatment vary con- 
siderably. Much of the Optics chapter is hardly more than a statement of 
physical principles, with little in the way of information about apparatus and its 
setting up. In the chapters on Heat, on Acoustics, and to a lesser extent in 
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those on Electronics and Atomic Physics, there is much more practical detail. 
Of intent, however, ‘ old, familiar experiments ’ are given little more than a men- 
tion and a reference, so as to make space for the less familiar and the new. 

The Manual has clearly been a labour of love by the friends and colleagues 
of Professor Ll. W. Taylor. In creating such a memorial, they must be con- 
gratulated on performing also a valuable service to Physics. It should serve as 
a guide and reference book to all who are involved in bringing up to date the 
advanced practical teaching of the subject, and in this instance ‘ advanced 
undergraduate ’ seems indeed to mean the same thing in both Britain and the 
United States. Vi BsC; 


CORRIGENDA 


Magnetoresistance of Copper, Silver and Gold, by M. G. PRrrestLEy, 1960, 
Phil. Mag., 5, 111. 


On p. 114 of this paper the last two lines of the table should read 
as follows: 
(110) 5 5 6-1 


Lae (211) 6 6 5-3 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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Relative temperature of platinum wire vs time when quenched into water. 
The temperature is measured by the resistance of the sample which is 
proportional to the voltage if the heating current is left on during the 
quench. 


L. H. GERMER et al. Phil. Mag. Ser. 8, Vol. 5, Pl. 26. 


Fig. 2 


C D 


Representative diffraction patterns from the (111) face of a clean nickel 
crystal (see table 1). (A) (533) reflections at 129 v : (B) (644) reflections 
at 200 v ; (C) (466) reflections at 262 v; (D) reflections in both of the 
principal azimuths at 120 vy. 


L. H. GERMER et al. PhilaMaguseracavole5, Pl. 27. 


= D 


Diffraction patterns representing a monolayer of nitrogen atoms on a (111) 
nickel surface, arranged in the 4 structure (see table 2). (A) 52 Vv, 
diffraction beams strong in B azimuth; (B) 57 Vv, diffraction beams 


strong in A azimuth ; (C) 65 v, diffraction beams strong in B azimuth ; 
(D) 77 v, diffraction beams in A and B azimuths equally strong. 


Cris GERMERVeual: Phil. Mag. Ser. 8, Vol. 5, Pl. 28. 


Fig. 5 


Diffraction pattern from unknown material upon crystal, before it had been 
thoroughly cleaned. Taken at 78 v, 1:39 A. 


Diffraction pattern from a monolayer of nitrogen atoms on a (111) nickel 
face, arranged in the 3 structure. 62 v (see table 2). 


L. H. GERMER et al. PhilsMagvisen cm Voleo, Ple29. 
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Diffraction patterns from nitrogen atoms arranged in the Q structure (see 
table 3). (A) 48v; (B) 75v; (C) 120v, pattern diffuse ; (D) 120Vv, 
pattern sharp after heating to 375°c. 


vice GERMER? et.al: Phil. Mag. Ser: 8) Volz5, PIa30: 
Fig. 10 


Diffraction pattern from a monolayer of CO molecules arranged in the 
4 structure. The pattern is as sharp as one would expect from a single 
two-dimensional crystal. 


A B 


Diffraction patterns from a heavy layer of oxygen on the crystal surface, so 
much oxygen that no patterns at all appear from the crystal itself. 
82v, 1:35 4. Before anneal even the oxygen gave no _ pattern. 
(A) after anneal at 200°; (B) after further anneal. Heating the 
crystal sufficiently to bring back the Laue beams to full intensity 
resulted in no further change of the oxygen pattern. Apparently the 
Q structure (see table 4). 


A. D. WHAPHAM and M. J. MAKIN Baile Magmoeras,.VOl.5, Pl.3l. 


Fig. 4 


Group of three cracks seen in a lithium fluoride specimen which had been 
irradiated with 4x 10 1 Mev electrons/em? and then compressed 0-7°, 
along [001]. The cracks which form on the usual (100) cleavage planes 
parallel to the compression axis, are initiated at slip-plane intersections. 
The slip planes are made visible by enhanced optical absorption produced 
by deformation of the irradiated crystal. (90 x .) 


Fig. 5 
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SLIP PLANES ING 


The upper crack of the group of three in fig. 4 at higher magnification. 
‘There is no sign of any component of the crack on any plane other than 


(100). (500 x.) 


A. D. WHAPHAM and M. J. MAKIN Phil) Mag. Sers 8, Volmo, Plaa2: 


Lithium fluoride single crystal annealed for + hr at 700°c after an irradiation 
dose of 2 x 10!® I Mev electrons/em?, showing the rectangular cavities 
produced by vacancy condensation. The diagonal lines lie on {110} 
slip planes and are thought to result from non homogeneous, radiation- 
induced expansion of the crystal lattice. 


Phils Magi set..8..VOleo; Pl. 33: 
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M. F. ASHBY and G. C. SMITH Phil. Mag. Ser. 8, Vol. 5, Pl. 34. 
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Fig. 3 
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